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SUMMARY 

The objective of the study was to develop and evaluate a locally 
homogeneous flow (LHFj model of spray evaporation which employs a second - 
order turbulence model. The LHF approximation is a simplified approach 
for modeling two-phase flows, where it is assumed that interphase trans- 
port rates are infinitely fast. The turbulence model was based on the 
Reynolds-averaged form of the conservation equations. The equations 
were solved using the GENMIX computer program. Solutions were obtained 
for both single and two-phase jets in stagnant air. 

Measurements to check the predictions were obtained for an isothermal 
air jet, an isothermal dense gas jet (sulfur hexafluoride), and an 
evaporating spray (Freon- 11, Sauter mean diameter 29 ym) . A single- 
channel LDA was used to measure mean velocities, turbulent fluctuations 
and.Reynolds stress. Composition was measured by isokinetic sampling 
and gas chromatography. Temperatures were measured with a fine wire 
thermocouple. Drop size distributions and liquid flux were measured 
by slide impaction. The evaluation also considered measurements reported 
in the literature. The evaluation proceeded systematically through the 
following flows: 

1. isothermal air jet in air 

2. heated air jet in air 

3. dense gas jet in air 

4. air jet in water 

5. evaporating spray in air 

The model provided reasonably good predictions of mean and turbulent 
quantities in the single phase flows, using constants in the turbulence 
model established during earlier work. The model provides a reasonably 
good estimation of radial variations of mean and turbulent quantities in 
the spray, but tends to overestimate the overall state of development 
in the axial direction. For example, in the fully developed region, 
predicted values of mean axial velocity and mixture fraction tended to 
be 10-20 and 40 percent below the measurements, respectively. 

Predicted drop life histories along the centerline of the jet were 
used to assess the effect of finite rate interphase transport on the 
LHF model predictions. The results suggest that a spray having a Sauter 
mean diameter less than 10 ym would be required for quantitative accuracy 
with the LHF model, for present test conditions. The need for such 
small drop sizes for quantitative accuracy with the LHF model is due 
to the rapid development of the flow resulting from the small injector 
diameter. Based on these findings, LHF models provide a useful qualitative 
picture of spray structure; however, more complex models allowing for 
finite interphase transport rates are required for quantitative accuracy 
in most practical sprays. 

The authors wish to acknowledge the assistance of H. Tamura, 
of the National Aerospace Laboratory, Kakuda Branch, Japan, While 
in residence at Penn State, Mr. Tamura contributed significantly to 
the development of the computer program used in this study. 



CHAPTER I 


INTRODUCTION 

1.1 General Statement of the Problem 

The combustion of liquid fuels accounts for a significant fraction of 
the current energy supply. For example, liquid petroleum products 
provided 46% of the energy supply in the U.S. in 1973, with a similar 
fraction being representative of world supplies (1). For combustion 
to occur, fuel and air must be intimately mixed. Almost without 
exception, mixing involves the use of an evaporating spray at some point 
in the combustion process. In some cases mixing is separated from 
combustion as in the carburetor of a spark ignition engine. However, 
in most cases, the evaporation and combustion of the spray occur 
simultaneously at the point of injection. The overall objective of the 
present study is to model the spray evaporation process with a simplified 
but widely used approach and evaluate the model with a series of well 
defined experiments. 

The spray evaporation process to be considered involves the vaporization 
of a well-atomized liquid jet. A spray is considered to be well- 
atomized if the liquid jet disintegrates into a mist of droplets 
immediately after injection and if the droplets are small in comparison 
to the injector diameter. For a constant area injector with no swirl, 
the ratio of the injector diameter to the diameter of the largest drops 
should be in the range of 10 to 100 to be considered well-atomized (2). 
After being formed at the injector exit, the droplets travel through 
the flow field until complete vaporization occurs. 

Figure 1 is a schematic representation of the spray evaporation 
process. The nozzle is a full -cone pressure atomized injector with 
no swirl. The surrounding environment is a stagnant inert gas which 
does not contain the evaporating substance as a vapor. 

The distribution of droplets in the spray is nonuniform with the smaller 
drops on the periphery and the larger drops near the centerline. The 
small drops readily exchange momentum with the gas causing the spray 
to entrain the surrounding gas similar to a gas jet. Near the injector 
exit, the velocity difference between the drops and the gas is the 
greatest, and the momentum of the liquid is transferred to the gas over 
an extended distance. Therefore, the development of the turbulent 
spray is retarded in comparison to a gas jet which affects the rate of 
entrainment and mixing characteristics of the flow. 

The behavior of the large drops is significantly different from that of 
smaller droplets. Small drops tend to follow the flow and are largely 
confined to fuel-rich eddies. The paths of larger droplets do not 
follow the turbulent eddies because of their greater inertia. The 
larger droplets pass through both injected fluid regions and regions composed 
largely of the entrained gas. The variation of composition and temperature 
between these regions significantly affects the evaporation rates of 
the individual droplets. 
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Figure 1 Schematic Diagram of the Spray Evaporation Process 
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The spray evaporation process is composed of several complex phenomena; 
the hydrodynamics of spray formation, the transport characteristics of 
the individual drops and the turbulent two-phase flow of the spray. 

All of these processes must be considered in a complete two-phase 
flow model. In order to model the spray the initial drop size and 
velocity distributions generated by the injector must be specified. 

This information is very difficult to acquire and is rarely available 
for practical sprays. The effect of turbulence on drop transport in 
actual spray conditions is not well understood and must be accounted 
for in a complete two-phase flow model. The effect of the droplets on 
the generation and dissipation of turbulence in a spray is not well 
defined. Solution of these problems will require extensive research 
into each area. 

One way to overcome many of the difficulties in a two-phase spray 
process is the use of a locally homogeneous flow model. The basic 
assumption of the locally homogeneous flow model is that the transport 
between the phases is infinitely fast. This implies that the droplets 
are in thermal equilibrium with the local gas environment (the gas 
immediately surrounding the droplet), and that there is no velocity 
difference (no slip) between the liquid and gas phases. 

The major advantages of the locally homogeneous flow model are as 
follows: 

1. The model requires little information concerning injector charac- 
teristics, e.g., the flow rate, the momentum, the composition and 
the temperature at the injector exit. 

2. The theoretical model of the flow is equivalent to that of a single- 
phase flow, bypassing the difficulties involved with correctly 
modeling interactions between the phases and minimizing the 
number of empirical parameters to be specified. 

3. Calculations with the model can employ the highly developed 
general purpose computer programs for single-phase flows, 
with only minor modifications. 

The major disadvantage involves inaccuracies due to failure of the 
assumption of infinitely fast transport rates between the phases. The 
difficulty becomes more important when drop sizes are large, when the 
rate of development of the process as a whole is fast (the process is 
completed near the injector) and when liquid densities are much greater 
than gas densities. 

A major objective of the present study will be to examine the capa- 
bilities of locally homogeneous flow models to predict spray evaporation 
and to develop a means of assessing their potential validity for a 
given application. 
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1.2 Previous Related Studies 


1.2.1 Locally Homogeneous Flow Models 

The locally homogeneous flow model has been applied to a wide range of 
two-phase systems. Both gas-liquid and liquid-gas processes have been 
modeled with the locally homogeneous assumption. These studies should 
provide the basis for a detailed model of spray evaporation. 

Weimer, Faeth and Olson (3) successfully modeled condensing gas jets 
injected into subcooled liquids using an integral model of the turbulent 
flow along with the locally homogeneous assumption. Vapor penetration 
lengths for condensing water, ethylene glycol and iso-octane jets were 
correlated irrespective of the vapor-liquid system. A variable density 
entrainment law was used to account for the mixing of the injected 
fluid and the surroundings . 

A turbulent air-water jet has been successfully modeled with the locally 
homogeneous assumption (4) . Measurements of the velocity, void 
fraction and momentum flux were completed throughout the jet. A 
locally homogeneous integral model was used to correlate the axial and 
radial distributions of the quantities. 

The locally homogeneous assumption has also been used to model a 
reacting gas-liquid system (5). This involved the combustion of a 
gaseous oxidizer jet in a liquid metal fuel. The locally homogeneous 
assumption was used to model the two-phase nature of the jet. A 
unified correlation of flame length, temperatures, and velocities was 
developed for the reacting two-phase system. The model also gave good 
results when applied to reacting single-phase jets and condensing gas jets. 

Newman and Brzustowski (6) applied the locally homogeneous assumption 
to a liquid carbon dioxide jet, with the liquid near the thermodynamic 
critical point, which was injected into a stagnant atmosphere of 
carbon dioxide and nitrogen. The locally homogeneous approximation is 
appealing near the critical region since the gas density and the liquid 
density are nearly equal. For this condition, the droplets quickly 
reach equilibrium with the gaseous environment. Mean velocities and 
temperatures were predicted quite well for the carbon dioxide and nitrogen 
system, however, the model was not tested against any other variable 
density systems. 

The locally homogeneous flow model has also been applied to a non- 
reacting Diesel fuel spray. Wakuri, et al, (7) conducted photographic 
measurements of fuel sprays with a high speed motion picture camera. 

A non-dimensional correlation was developed for the spray length as 
a function of the density ratio of the liquid to the gaseous environment, 
inlet velocity, initial diameter, and time. It was concluded that for 
their experimental conditions the relative velocity difference between 
the phases was negligible and therefore the two phases could be treated 
as a homogeneous mixture. 

There have been several attempts to extend the locally homogeneous 
assumption to the combusting spray in Diesel engines (8-12). However, 
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due to the transient three-dimensional nature of Diesel combustion, 
convincing evaluation of the models has not been achieved. The 
correlations developed by these studies are limited in their applica- 
tions to the design of Diesel engines and have not been extended to 
consider other spray combustion processes. 

Thring and Newby (13) analyzed the spray combustion process in a 
reacting oil spray. The two-phase spray was taken to be equivalent to 
a single-phase jet of the same momentum and equal stoichiometry. 

Qualitative agreement was obtained for the prediction of overall 
quantities, such as the spray length, based on the locally homogeneous 
assumption. 

The results of several other experimental investigations indicate 
that the locally homogeneous model can, at times, be valid for reacting 
two-phase sprays. Onuma and Ogasawara (14, 15) and Chigier and co- 
workers (16, 17) have investigated several types of two-phase reacting 
sprays at atmospheric conditions. In these studies the spray structure 
was defined by the gas temperatures, droplet sizes, velocities, and 
concentrations. It was concluded that droplets in a flame do not burn 
individually but that the fuel vapor diffuses from a region of droplet 
evaporation and reacts in a way similar to a single-phase reacting jet. 

However, the results of more recent investigations indicate that the 
locally homogeneous assumption cannot be successfully applied to all 
two-phase reacting sprays. Khalil and Whitelaw (18) measured the velo- 
city, turbulence intensity, temperature and droplet concentration in a 
hollow-cone spray at atmospheric pressure. The flame length for a 
spray with a Sauter mean diameter of 45 uni was larger than the flame 
length of the simulated gaseous diffusion flame. The difference indicates 
the assumption of a locally homogeneous flow, which does not provide for 
dropsize effects, can be invalid for sprays having a Sauter mean dia- 
meter as small as 45 ym. Recent calculations by Khalil (19) show better 
agreement between theory and experiment; however, this was obtained 
for selective cases by introducing several empirical parameters which 
have not been systematically evaluated. 

Failure of the locally homogeneous model has been illustrated for pressure- 
atomized reacting sprays even at elevated pressures. Shearer (20) 
compared experimental measurements of flame shape and spray length with 
theoretical predictions for a full cone spray at ambient pressures of 
100 to 9000 kPa. A locally homogeneous model, in conjunction with an 
integral model for the flow, correctly predicted the trends of the data, 
but underestimated the extent of the spray and flame boundaries by 30 
to 50 percent. Although the Sauter mean diameter of the spray was 
approximately 30 pm, it was concluded that there was significant slip 
velocity between the liquid and gas. 

In summary, the locally homogeneous flow model is based on the assumption 
that all phases have the same velocity and are in thermodynamic 
equilibrium at each point in the flow. This requirement is most easily 
met when the dispersed phase is a gas and the continuous phase is a 
liquid due to the relatively low inertia and thermal capacity of the 
gas phase. Therefore, locally homogeneous flow models have been reasonably 
successful for describing processes of condensation and reaction when gas 
jets are injected into liquids. 
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In contrast, the ability of the locally homogeneous model to predict 
the behavior of an evaporating or combusting spray has not been clearly 
defined. Qualitative agreement between predicted results and the experi- 
mental results has been obtained, however, the locally homogeneous model 
tends to overestimate the rate of development of the flow. An in-depth 
evaluation of the locally homogeneous assumption has not been conducted. 

1.2.2 Two-Phase Flow Models 

A number of two-phase flow models have been developed to predict the 
behavior of evaporating or combusting sprays. All of the two-phase 
models attempt to treat the interaction between phases; however, the 
models vary widely in regard to the thoroughness with which the inter- 
actions are treated. 

A major distinguishing feature of the present two-phase flow models is 
the consideration of the relative velocity (slip) between the phases. 

One group of models (21-24) ignores slip between the phases as opposed 
to more complete models (25-32) where slip is considered. If slip is 
ignored, the models are similar to the locally homogeneous flow models, 
except that finite heat and mass transfer rates are considered between 
the phases. Ignoring slip considerably simplifies the calculations; 
however, such models involve limitations in accuracy similar to locally 
homogeneous flow models. 

In most two-phase flow models it is assumed that the spray is dilute. 

This implies that single drop models (33-42) can be used to describe 
drop transport rates, that the presence of the drops has a negligible 
effect on the constitutive equations describing the transport character- 
istics of the gas phase, and that drop collisions can be ignored. 

The void fraction requirement of all of these assumptions is not well 
defined. For monodisperse particles, it is recognized that single 
particle transport formulas apply reasonably well if the center-to-center 
distance between two particles is greater than two particle diameters (43). 
If the particles have equal spacing, this requirement implies a void 
fraction greater than 0.92. 

For a reacting spray with no swirl in a stagnant environment, void 
fractions in the range of 0.92 are reached approximately 6 to 12 injector 
diameters from the injector (20). This suggests that the dilute spray 
assumption is adequate for the bulk of the two -phase region. However, 
the near injector region requires a different formulation. 

The manner in which the initial conditions of the spray are specified 
has a strong bearing on the need to treat the low void fraction region. 
Two-phase flow models which consider slip require size and velocity 
distributions as a function of position at some point in the flow. Few 
measurements of these quantities have been conducted due to current 
experimental limitations. 

Model development, however, has far outstripped systematic efforts at 
model validation. Where predictions have been compared with experiments, 
injector characteristics have not been defined adequately for a two-phase 
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flow model and measurements within the flow have been limited. The 
models are very complex, involving numerous assumptions and empirical 
parameters. Convincing validation will require careful testing of 
injectors and simple well defined experiments. Present two-phase flow 
models are not capable of accurate, a priori , predictions of evaporating 
or combusting sprays (2) . 

1.2,3 Turbulence Models 

Second-order turbulence models have been applied to a wide variety of 
single-phase non-reacting flows. Launder, et al., (44) developed a 
turbulence model based upon the turbulent kinetic energy and the dissi- 
pation rate. Solution of the turbulence model was accomplished by intro- 
ducing the equations for the turbulent kinetic energy and the dissipation 
rate into the GENMIX computer program developed by Patankar and Spalding 
(45) . The predicted values for the mean velocity and the turbulent 
kinetic energy were in good agreement with most of the experimental 
results. However, it was suggested that the analytical model would not 
adequately represent flows with large density gradients. 

A second-order turbulence model has also been successfully applied to 
single-phase reacting jets by Lockwood and Naguib (46). Turbulence was 
modeled using differential equations for the turbulent kinetic energy, 
the dissipation rate, and concentration fluctuations. A Reynolds- 
averaging procedure which neglects the effect of density fluctuations 
was used to model the conservation equations. A clipped Gaussian 
probability density function was employed to represent the relationship 
between mixture fraction and other properties. 

Lockwood and Naguib (46) used a systematic approach to evaluate the 
model. The predictions were compared with data for isothermal inert 
flows, non-isothermal inert flows and reacting flows, Good agreement 
was obtained between predictions and experimental results for all three 
cases, however additional data is required to completely evaluate the 
model. In particular, more information is needed for highly non-isothermal 
inert flows in order to fully evaluate the effect of large density 
differences in the flow. 

A second-order turbulence model which includes the effect of density 
fluctuations has been successfully applied to single-phase reacting 
jets by Gosman, Lockwood and Syed (47). Horizontal and vertical diffusion 
flames were modeled and the theoretical predictions were compared with 
measurements. The model employed equations for the conservation of mass, 
momentum, energy and mixture fraction in conjunction with a second-order 
turbulence closure. The turbulence was modeled with equations for the 
turbulent kinetic energy, dissipation rate, and concentration fluctua- 
tions. Terms involving the density fluctuations were retained in portions 
of the model, although complete variable density modeling was not achieved. 
The equations for three dimensional flow were solved using the TEACH-3P 
computer program (48) . 

The Gosman, Lockwood and Syed (47) model has only been tested against 
measurements of temperature and oxygen concentration from a single 
experiment. The results for the vertical reacting jet compare moderately 
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well with the predictions, however, additional evaluation of the model 
is required before any definite conclusions concerning the validity of 
all aspects of the model can be reached. 

The combination of a second-order turbulence model and a Favre-averaging 
procedure has been successfully applied to a single-phase reacting 
jet. A Favre-averaged quantity represents a density weighted average. 

In this approach, the fluctuating density effects are included in the 
density weighted average. Kent and Bilger (49) used Favre-averaged 
equations for the conservation of mass, momentum, and mixture fraction 
and a turbulence model based on the turbulent kinetic energy, dissi- 
pation rate, and concentration fluctuations to model a reacting single- 
phase jet. With their turbulence model in the Favre-averaged form, the 
equations to be solved are identical to the incompressible flow equations. 

The Favre-averaged model (49) was tested against measurements of 
velocity, temperature, and pollutant concentrations for a reacting jet 
in a co-flowing air stream. The experimental measurements and theoreti- 
cal predictions compare very favorably. Evaluation of the model, 
however, was limited to examining a single data set. 

The effectiveness of several turbulence models has been examined in a 
variable density planar shear layer by Parker and Sirignano (50) . 

The study included a mixing length model, a Favre-averaged second order 
turbulence model, and a conventional Reynolds-averaged second order 
turbulence model which ignores density fluctuations. For the case of a 
variable density planar shear flow the Favre-averaged model provided a 
better prediction of the boundary layer than the Reynolds-averaged 
model. However, the mixing length model provided a better prediction 
of the flow than either of the second-order turbulence models. 

Second-order turbulence models have been applied to a wide variety of 
single-phase flows. However, only the Lockwood and Naguib (46) model 
has been systematically studied for a wide range of variable density 
flows. There is also evidence to suggest that the effect of density 
fluctuations is significant and should be considered. 

1.3 Specific Objectives of the Study 

The availability of an accurate two-phase flow model for sprays would 
provide a most useful design tool. Models of this type, however, require 
systematic validation, which is difficult at the present time due to 
experimental problems in accurately describing initial drop size and 
velocity distributions which are needed in the analysis. Locally 
homogeneous flow models provide a simpler approach, where the need for 
detailed injector characteristics is eliminated. While the locally 
homogeneous flow method provides a logical first-step toward the goal 
of validated two-phase flow models, the previous discussion indicates 
that systematic validation of even this simplified approach has not yet 
been reported. 
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Based on this assessment of current knowledge, the objective of the 
present investigation is to develop and evaluate a locally homogeneous 
model of spray evaporation. Following the recent studies, a second- 
order turbulence model will be employed in the analysis. The turbulence 
model used by Gosman, Lockwood and Syed (47j has been subjected to the 
most systematic evaluation of the various methods reported to date; 
therefore, this procedure will be employed in the present study. Aside 
from comparing the predictions of the model with spray experiments, effects 
of slip will also be evaluated directly by completing drop life-history 
calculations using the predictions of the model to supply the local 
conditions surrounding the drop. 

The evaluation of the model will consider a wide range of conditions for 
both single and two-phase flows. This will include single-phase constant 
density jets, single-phase variable density jets, two-phase gas-liquid 
jets and evaporating sprays. In this manner, the various parameters in 
the model can be established in single phase flows, prior to evaluation 
of the locally homogeneous flow approximation in the two-phase flows. 

The single-phase constant density jet is the least complicated flow to 
model. In this case, the predictions can be compared with existing air 
jet measurements reported by Wygnanski and Fiedler (51), Becker, et al., 
(52) and Hetsroni and Sokolov (53). Experiments in this flow will also 
be conducted in this investigation in order to establish the validity 
of present experimental techniques by comparison with the earlier 
measurements. 

Data for single-phase variable density jets are more limited (54). 
Therefore, additional measurements in this flow will be undertaken in 
the present investigation. This will involve a sulfur hexafluoride gas 
jet, which has a high density in comparison to air due to its higher 
molecular weight, injected into stagnant air. 

Tross (4) provides data on velocity and composition profiles in the 
flow resulting when an air jet is injected into a stagnant water bath. 

These results will be employed to test the locally homogeneous flow 
model under conditions where this assumption is generally satisfied. 

Data for an evaporating spray will be generated in the present investiga- 
tion. The configuration will consist of a Freon- 11 spray produced by 
an air atomizing injector without any swirl. The flow will be injected 
into stagnant air. 

The specific objectives of the present study can be summarized as follows: 

1. Develop a model of the spray evaporation based on the locally homo- 
geneous assumption in conjunction with a second-order turbulence model. 

2. Experimentally measure mean and fluctuating velocities (u,v,u’,v-, 
w',u’v') throughout single-phase constant density and single-phase 
variable density jets as well as a two-phase evaporating spray. 

3. Experimentally measure the injected fluid concentration throughout 
single-phase variable density jet and two-phase evaporating spray. 
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4. Experimentally measure the gas temperature and drop size throughout 
an evaporating spray. 

5. Compare the theoretical predictions of the model with the experi- 
mental measurements and previous studies. 

6. Develop a model for the behavior of a single drop in the evaporating 
spray to determine the effect of slip in the evaporating spray. 


CHAPTER II 

EXPERIMENTAL APPARATUS AND PROCEDURE 

2 . 1 Introduction 

The experimental apparatus was used for measurements in both single 
and two-phase evaporating sprays. The flows examined during the 
present study are as follows: 

1. Isothermal air jet in air 

2. Isothermal dense gas jet in air (sulfur hexafluoride) 

3. Evaporating spray in air (Freon-11). 

In all cases, the jet issued vertically downward into stagnant air. 

A variety of measurements were performed to define the flow. Jet flow 
rates and the total momentum of the jet were determined to define 
initial conditions. Profiles of mean gas velocity, velocity fluctuations, 
and Reynolds stress were measured for all flows. Mean concentration 
measurements were completed for the isothermal dense gas jet and the 
evaporating spray. Mean temperature and dropsize distribution profiles 
were measured for the evaporating spray. 

2 . 2 Experimental Apparatus 
2.2.1 Test Apparatus 

The nozzle was positioned near the top of the test stand and oriented 
vertically downward. The test stand was an area 2.0 m high by 0.92 
m square enclosed with a 16 mesh screen. Two traversing mechanisms 
were used to position the nozzle in a horizontal plane. 

An air-atomizing nozzle was used throughout the investigation. The 
nozzle was manufactured by Spraying Systems Company (Model 1/4 J 2050 
fluid nozzle and 67147 air nozzle). The outlet diameter of the nozzle 
was 1.194 mm. 

The fluid supply system consisted of a gas supply portion and a liquid 
supply section. A schematic of the fluid supply system is shown in 
Figure 2. The same fluid supply system was used for all three flows. 

The gas supply system controlled the air flow rate for the evaporating 
spray and the isothermal air jet and the gas flow rate for the variable 
density single-phase jet. The gas flow rates were metered and controlled 



20 MPa 
AIR SUPPLY 



Figure 2 Spray Evaporation Test Apparatus 
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with a critical-flow orifice and pressure regulator combination. The 
pressure regulator was a Matheson model 4 regulator with a 0-20.68 MPa 
output capacity. The pressure upstream of the critical-flow orifice 
was measured with a Heise absolute pressure gauge with a 0-5,17 MPa 
capacity. The orifice was calibrated with a Precision Scientific 
Company wet-test meter (0.283 £/rev.]. 

The liquid supply system was only required for the evaporating spray 
tests. The liquid was contained in a feed tank with a calibrated 
sight glass. The liquid flow rate was measured by timing the fall of 
the liquid level in the sight glass. The fluid was delivered to the 
nozzle by pressurizing the feed tank with air. A Matheson model 8 
regulator was used to control the pressure. The pressure was measured 
with a 0-0.41 MPa gauge manufactured by Ashcroft Corporation. A 
liquid supply tank was connected to the system to supply additional 
liquid. 

2.2.2 Velocity Measurement 

Velocity measurements in the single-phase jets and the two-phase 
evaporating spray were conducted with a single-channel laser Doppler 
anemometer (LDA) . An equipment list for the LDA is supplied in Table 1. 
The sending and receiving optics had a focal length of 595 mm with a 
4.7® angle between beams. The aperture diameter of the photodetector 
was 0.256 mm. The signal was focused on the photodetector with a 200 mm 
lens. This produced an ellipsoidal measuring volume 6.224 mm in length 
with a diameter of .763 mm. A schematic of the LDA arrangement is shown 
in Figure 3. 

Oil drops were used as the seeding material in the flow. The exhaust 
from a vacuum pump operating at its maximum flow rate contained a large 
number of suspended oil particles. These oil droplets were used to seed 
the surrounding gas. The average diameter of the particles was 0.618 
ym at a concentration of 2,8 x lO^O particles/m^. The injected flow 
was not seeded, which limited measurements to the region far from the 
injector where a large number of oil particles had been entrained. 

Addition of the seeding particles changed the density and specific heat 
of the air less than 0.01 percent. 

The seeding characteris tics of the air was analyzed with a Royco model 
230 light scattering particle counter. A 100 to 1 sample diluter was 
used to bring the concentration of the air sample within the range of 
the instrument. The samples were examined for one minute periods at 
0.283 liters per minute flow rate. 

Mean and fluctuating velocity components were measured in the axial, 
radial and tangential directions at several axial distances from the 
injector for all three test conditions. The measurements were repeatable 
within 10 percent for different traverse^. Bec ause the LDA system was 
limited to a single channel, u and u'^, v and v’ ^ , and w' ^ were 
measured during three separate traverses of the flow with a different 
beam orientation for each traverse. An integration period o f one minute 
was used to determine the mean quantities. A measurement of u* v ’ may also 
be derived by combining the results of several traverses (55) . The 
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Table 1 


LDA Equipment List 


Component 

Manufacturer 

Model 

Helium-Neon Laser 

Spectra Physics 

125A 

Integrated Optics 

Thermo-Systems 

900 

Frequency Shifter 

Thermo-Sy s t ems 

980 

Photodetector 

Thermo- Sys terns 

960 

Tracker 

Thermo-Sy s terns 

1090 

RMS Voltmeter 

Thermo-Sy s terns 

1076 

Dual Beam Oscilloscope 

Tektronix 

912 


Integrating Voltmeter 


Hewlett-Packard 


240IC 





INJECTOR 




Figure 3 Schematic Diagram of the LDA System 
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evaluation of the turbulence quantities from these measurements is 
described in Appendix A. 

2.2.3 Concentration Measurement 

Mean concentration measurements of the injected fluid were conducted 
for the variable density single-phase jets and the two-phase evaporating 
spray. Mean concentration measurements across the entire radius of the 
flow were obtained at axial positions of 203, 406 and 609 injector 
diameter from the injector. 

The sampling probe tip is illustrated in Figure 4. The flows were 
isokinetically sampled by applying suction to the sampling probe and 
measuring the sampling flow rate. The sampling flow rate was adjusted 
to match the local flow rate determined from the velocity measurements. 

The sampling flow rate was measured with a Precision Scientific wet- 
test meter (0.0283 £/rev.). The tip of the sampling probe could be 
heated to prevent water in the air from condensing and freezing at the 
probe inlet during evaporating spray tests. 

The samples were analyzed with a Perkin Elmer Model 880 gas chromato- 
graph using a hot-wire detector. A 5 ml gas tight syringe was used 
to transfer the sample from the sampling port to the gas chromatograph. 

The separating column for sulfur hexafluoride was stainless steel, 

6.35 mm O.D. by 122 cm in length packed with 80-100 mesh silica gel. 

The column was placed in an ice bath and operated at 0° C. The detector 
was operated at 250° C. Helium was used as the carrier gas and maintained 
at a flow rate of 15 cc/min. 

The chromatograph was calibrated for sulfur hexafluoride-air mixtures 
using a commercially prepared sample purchased from Matheson Gas 
Products. The calibration procedure indicated that in the range of 
interest the peak area was a linear function of concentration. 

The separating column for Freon-11 was stainless steel, 3.175 mm O.D. 
by 183 cm in length packed with PORAPAK Q (80-100 mesh) . This column 
was placed in the chromatograph oven and maintained at 150^ C. The 
detector was maintained at 250° C. Helium was used as the carrier gas 
at a flow rate of 20 cc/min. Various mixtures of Freon-11 and air and 
a commercially prepared sample of Freon-11 and helium (purchased from 
Matheson Gas Products) were used to calibrate the chromatograph. The 
calibration procedure indicated that the peaks were symmetrical and the 
concentration was a linear function of peak height for the concentration 
range of interest. 

2.2.4 Temperature Measurement 

Gas temperatures were measured at several axial locations in the spray 
with a fine wire thermocouple. Temperature measurements were conducted 
across the entire radius of the spray at several axial positions. The 
temperature measurements can only be considered qualitative estimates of 
the gas temperature because the effect of droplet impingement on the 
thermocouple was not considered. 



CERAMIC 

INSULATION 

HEATER 


GLASS TUBE 


Figure 4 Sketch of the Sampling Probe 
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An illustration of the thermocouple probe is shown in Figure 5. The 
probe was manufactured from 0.076 mm diameter Chrome 1-Alumel wire and 
was spot welded onto 0.75 mm diameter wire leads of the same material. 

The reference junction was placed in a reference cell manufactured by 
Hy-Cal Engineering maintained at 0° C. A Hewlett-Packard integrating 
digital voltmeter, model 240 IC, was used to integrate the signal over 
a one minute period, in order to determine the mean temperature. 

2.2.5 Drop Size and Liquid Flux Measurements 

An inertial impaction method was used to determine the droplet size 
distribution and the liquid mass-flux. Small glass slides were coated 
with magnesium oxide and exposed to the two-phase spray for a specified 
time period. As each drop hits the slide, it leaves an impression in 
the oxide coating, proportional to its size. The slide was examined 
under a microscope and each drop impression was sized and counted. This 
permitted determination of the dropsize distribution over a cross-section 
of the spray and the point-to-point variation in the dropsize across 
a spray cross-section. A measurement of the liquid mass-flux was also 
determined. These measurements were performed at several axial locations 
to indicate the axial variation in dropsize distribution and the liquid 
mass-flux. 

A shutter mechanism was designed to expose the slide to the spray for 
a given time period. A sketch of the shutter mechanism is shown in 
Figure 6. Compressed air forces the shutter to move across the slide. 

The shutter speed was timed at 20 milliseconds. The collection effi- 
ciency of the shutter mechanism may be calculated for a range of droplet 
sizes for a given relative velocity (56). This was used to correct the 
number of droplets collected. 

The droplet size distribution was measured at several axial locations. 

At each axial location, the droplet distribution was measured across 
the entire width of the spray. The Sauter mean diameter and the average 
droplet diameter are determined directly from the droplet size distri- 
bution. The liquid mass flux was determined by dividing the total 
liquid mass collected by the^ collection area and the exposure time. 

2.2.6 Je t Momentum Measurement 

The momentum of each flow was determined by measuring the axial force on 
an impingement plate held near the exit of the nozzle (100 mm square 
plate, 25 mm from the nozzle). A Unimeasure 80 force transducer was 
used to measure the force. The transducer was calibrated, in turn, 
by placing known weights on the plate. The initial velocity was then 
determined from the jet momentum and the mass flow rate. 

2.3 Experimental Conditions 

The air for the system was supplied by a reciprocating type Ingersoll- 
Rand compressor. The compressed air was filtered with a Matheson type 
451 filter to remove oil and particles greater than 5 pm. 
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Figure 5 Sketch of the Thermocouple Probe 



2.54 


SUPPORT ROD 6.35 DIA. 


GLASS SLIDE 






DIRECTION 
OF TRAVEL 





GUIDE ROD 
6.35 O.D. TUBE 


-^TO AIR 
SUPPLY 


Figure 6 Sketch of the Droplet Impactor 




20 


The sulfur hexafluoride was produced by Matheson Gas Products. The 
gas was contained in a cylinder under its own vapor pressure (2200 
kPa at 294 K) . The gas was certified purity grade with a listed purity 
of 99.8 percent. 

The Freon-11 was supplied by E.I. DuPont de Nemours and Company. The 
purity of the liquid was 99.9 percent. 

The test conditions for the three cases are summarized in Table 2. 


CHAPTER III 

THEORETICAL CONSIDERATIONS 
3 . 1 Description of Model 

The theoretical model considers a steady, a xi symmetric, turbulent jet 
in an infinite, stagnant media. The overall configuration of the flow 
is illustrated in Figure 1. The jets to be considered are as follows: 

1. Isothermal air jet in air. 

2. Heated air jet in air. 

3. Dense gas jet in air. 

4. Air jet in water. 

5. Evaporating spray in air. 

The analysis generally follows the procedure developed by Lockwood and 
coworkers (46, 47), since this approach has been successful in earlier 
analyses of single-phase jets. The contribution of the present analysis 
is consideration of two-phase jets using the locally homogeneous flow 
approximation . 

The locally homogeneous flow approximation for a gas-liquid jet or 
evaporating spray implies that liquid and gas velocities are the same 
and that the two-phase mixture is in thermodynamic equilibrium at each 
point in the flow. This means that the temperature of both phases is the 
same and that phase equilibrium is maintained, i.e., the chemical 
potential of each species is the same in both phases. 

In order to satisfy the locally homogeneous flow approximation, the rates 
of transport between the phases must be fast in comparison to the rate 
of development of the flow as a whole. The validity of the approximation 
can best be assessed when some information is available concerning 
conditions within the flow. In the present study, the locally homogeneous 
flow model was used as a first estimate of the flow, followed by particle 
trajectory calculations in order to evaluate transport rates between the 
phases . 

Utilizing the assumption of locally homogeneous flow, the remainder of 
the analysis parallels models for single-phase flow. Typical of 
analyses of low Mach number jets, it is assumed that the boundary layer 
approximations apply, and that viscous dissipation effects and kinetic 
energy are negligible. Due to the relatively low temperature levels and 
small temperature differences encountered in the measurements to be 
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Table 2 



Summary 

of Test Conditions 



Case 

1 

2 

3 

Injector Fluid 

Air 

Sulfur Hexafluoride 

Freon-11 

Injector Flow Rates 
(g/s) 




Gas 

0.467 

0.103 

0.225 

Liquid 

0 

0 

1.548 

Injector Pressures 
(kPa) 




Gas 

203 

186 

211 

Liquid 

— 

— 

138 

Jet Momentum 
(N) 

0.100 

0.0093 

0.132 

Sauter Mean Diameter^ 
(ym) 

— 

— 

29 

Initial Velocity 
(m/s) 

214.1 

90.30 

74.45 


= 296 K, P =97 kPa 

00 ^ CJO 

^Measured 203 mm from injector. 
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considered with the analysis, radiation is also ignored. 

Following Lockwood and coworkers (46, 47) among others, it is assumed 
that the exchange coefficients for all species and heat are the same. 

This implies equality of both the laminar and turbulent components of the 
exchange coefficients. This is a reasonable approximation for single- 
phase flows (46, 47, 57, 58), but is more questionable for gas-liquid 
mixtures. For example, even at the locally homogeneous flow limit, the 
Iciminar diffusivity of small drops is much smaller than gas molecules, 
although turbulent diffusion rates are nearly the same (59) . However, the 
approximation is reasonable, since the effect of laminar diffusion in 
turbulent flow is not thought to be large, at least for mean properties 
(58). 

The assumption of equal exchange coefficients for all species and heat 
implies that, f, the mixture fraction (defined as the fraction of mass 
at a given point which originated from the injector) is a passive scalar 
or conserved property of the flow. In combination with the other assump- 
tions discussed earlier, this implies that the properties at each 
point in the flow correspond to the thermodynamic state attained when 
an amount f of injector fluid and (1-f) of ambient fluid, at their 
initial states, are adiabatically mixed at the ambient pressure of the 
j et. 

The turbulent flow model is based on the approach used by Gosman, Lockwood 
and Syed (47) for single-phase jets. This involves the solution of 
Reynolds-averaged conservation equations. The transport of mean quantities 
is given by conservation equations for mass, momentum and mixture fraction 
(equilibrium thermodynamics gives all other properties, temperature, 
density and composition, as a function of mixture fraction) . Turbulence 
characteristics are based on a second-order model requiring the solution 
of model equations for turbulent kinetic energy, dissipation and con- 
centration fluctuations (k-e-g model) . While buoyancy is considered in 
the mean equations, its effect on turbulence production and dissipation 
is ignored. This general procedure has been successful for a variety of 
turbulent single-phase jets (46, 47, 57, 58, 60-62). 

3. 2 Equation of State 

The relationship between mixture enthalpy, composition, temperature 
and density and the mixture fraction is provided by the equation of 
state. Each type of flow to be considered requires a separate 
equation of state. 

The experimental results to be compared with the analysis were all 
obtained at relatively low pressures. Therefore, all gases are assumed 
to be ideal gases. Since combustion is not treated, chemical reactions 
are ignored. Finally, it should be recalled that the mixing process 
is adiabatic and occurs at constant pressure, under the assumptions of 
the present analysis . 

Considering N species in the flow, the expression for the composition 
of the mixture is: 
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(3.1) 


where each species may exist in both the gaseous or liquid state 


Y. = Y-. + Y . . 

1 fi gi 


The enthalpy of the mixture can be expressed as 


h = h^f = h„(i-f) 


(3.2) 


(3.3) 


where at any condition 


h = I (Y h + Y .h .) i = 1, . . . ,N 


(3.4) 


The density of the mixture is given by 


P = I Y ,V .) 

^ fi fi gi gi^ 


-1 


i = 1,.. .,N 


(3.5) 


where and v . are the partial specific volumes of species i in the 
liquid and gas phases. 

Given the relationships between enthalpy and density of each species 
and the composition, temperature and pressure. Equations (3.1) - (3.5) 
are sufficient to describe the composition, temperature and density 
of the mixture. The relative composition of the gas and liquid phases 
is obtained from Equation (3.2) and the requirement that the chemical 
potential of each species must be the same in both phases. This 
requires that each type of flow must be examined separately. 

A gas jet injected into a stagnant gas environment is the least compli- 
cated flow to be examined. The present study considered three gas-gas 
systems: an isothermal air-air jet, a heated air-air jet, and a sulfur 
hexafluoride-air jet. All gases were modeled as ideal gases with a 
constant specific heat. Details of the gas-gas equations of state are 
presented in Appendix B.l - B.3. 

An equation of state for the two-phase air-water system investigated 
by Tross (4) has been formulated in the present study. The flow was 
considered to be isothermal and the air was assumed to behave as an 
ideal gas. It was assumed that the effect of water vapor was negligible 
and that no air was dissolved in the water. The details of the air-water 
equation of state are described in Appendix B.4. The equation of state 
for the air-water system is shown in Figure 7. The mass fractions of the 
air and water vary linearly with the mixture fraction. However, the 
variation of density with mixture fraction is very non-linear. 

The most complicated flow considered was the case of an evaporating 
spray. It was assumed that all gases behaved as ideal gases and that no 
air was dissolved in the liquid phase. The equation of state for the 


MASS FRACTION yodcg/m^) 
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Figure 7 Equation of State for the Air-Water System [4] 
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evaporating spray is shovm in Figure 8. In this case, both air and Freon 
11-liquid and vapor leave the injector, while the surrounding gas con- 
tains no Freon-11. The total mass fractions of air and Freon-11 are 
linear. The presence of liquid, however, causes the non-linear 
behavior of mixture temperature and density. Details of the evaporating 
spray equation of state are described in Appendix B.5. The physical 
properties used for all the flows are summarized in Appendix B.6. 

3 . 3 Probability Density Function Model 

The mean value of any scalar quantity 6 (other than f, k, e and gj can 
be determined from its variation with f if the probability density 
function of f, P(f) is known as a function of position in the flow. 

Given P(f), the mean value of any scalar is given by (46) 


0 


1 

0(f) P(f) df 

0 


(3.6) 


Following Lockwood and Naguib (46), a clipped Gaussian probability 
density function was assumed in the present study. This distribution is 
represented by a Gaussian function for the range 0 ^ f ^1, with the 
tails of the distribution being represented by delta functions at f = 0 
and 1, respectively. This distribution is represented by the following 
equations 


“ T72— C- -r 0<f<l (3.7) 


a(2lT) 


2 ^ a 


P(0) = A = 


a(2TT) 


T/2- t 4“ 


(3.8) 


P(l) = B = 


a(2TT) 


1 p 1 . f-y .2-, 

-TTr C- — (-1)—) J 


2 ^ a 


(3.9) 


and is illustrated in Figure 9. 

The most probable value of the distribution, y, and the variance, a, 
can be determined by noting that 


f = 


g = 


f P(f) df 


(f-f)^ P(f) df 


(3.10) 


(3.11) 




f 


Figure 8 Equation of State for the Evaporating Freon-11 Spray 







Figure 9 Probability Density Function for Mixture Fraction 
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Both f and g are known from the integration of the conservation equations. 
Therefore, Equations (3.10) and (3.11) provide two implicit equations 
to solve for a and y. This completes the specification of the probability 
density function. Values of p, T, Y. , etc., can then be obtained by 
integrating Equation (3.6), where 0 (^) is specified by the equation of 
state. 


3.4 Governing Equations 

Employing the assumptions discussed in Section 3.1, the Reynolds -averaged 
form of the conservation equations of mass, momentum, mixture fraction, 
turbulent kinetic energy, turbulent dissipation rate and concentration 
fluctuations are obtained as described in Appendix D. The final form of 
the equations are as follows: 

Conservation of Mass 


9pu 1 d 

9x r 8r 


(rp V ^) = 0 


(3.12) 


Conservation of Momentum 
D(u) = a(p^ - p) 


(3.13) 


Conservation of Mixture Fraction 
D(f) = 0 


(3.14) 


Conservation of Turbulent Kinetic Energy 


D(k) = y^ (^^) - pe 


(3.15) 


Conservation of Turbulent Dissipation 


D(e) = C 


C— 
^ 9r ^ 



(3.16) 


Conservation of Concentration Fluctuations 


D(g) 


C 


g 



^f_ 2 
9r '' 



where 


0 

P 


= p V + p' v' 


(3.17) 


(3.18) 
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and 


80 0 90 

D(4>) = P U + p V 


for 0 = u, £, k, e or g. 


3r 


Cr 


90 

"Ir" 


(3.19) 


The boundary conditions for these equations are: 


0=0 (3.20) 

f = 1 , g = 0 

(3.21) 

The parameters Kq and Gq are assumed to have the values appropriate for 
fully developed flow in tubes. 

The turbulent viscosity is obtained from the following constitutive 
equation : 

Pt = Pk^/£ (3.22) 

Equations (3.12) - (3.22) represent the final set of equations to be 
solved in this analysis. This set is applicable to all the flows 
considered in the present study. 


r = 0, 


90 

9r 


= 0 


r oo 


The initial conditions of the flow are 


X = 0 , r < d/2 , u = M^/m^ , 


k = kg . e = Cg 


3. 5 Numerical Solution 

3.5.1 Generalized Computer Code 

A general purpose computer code, GENMIX, developed by Patankar and 
Spalding (45) for boundary layer flows, was used to solve the governing 
equations. In order to solve the governing equations using GENMIX (45), 
they must be transformed into a coordinate system based on the axial 
distance x and a dimensionless stream function w. The dimensionless 
stream function is defined as: 


LO 


0-0 


(3.23) 


where 0, and 0^ are the values of the stream function at the inner and 
external boundaries of the flow. 



The stream function is defined as: 


8x 


and 

9r 


0 

-r p V 


r p u 


(3.24) 


(3.25) 


The use of the stream function insures the solution of Equation (3.12). 

At the inner and external boundaries the stream function can be expressed 
as 


9x 




(3.26) 


and 

9x 


-r^m 


EE 


(3.27) 


where m " and m^" are the mass transfer rates across the inner and 
external boundaries. 


These transformations permit Equations (3.13) through (3.18) to be 
put in the general form 


MA B(.) 


9o) 


9cjo 


(C 


9w ^ 


+ D 


(3.28) 


where (j) represents the dependent variables u, f, k, e and g. The other 
terms in Equation (3.28) are defined in Table 3. This set of equations 
is integrated from the injector to a value of x equal to 550 diameters 
downstream for values of u) ranging from zero to one. 

The GENMIX program (45) has the capability for as many as 30 cross- 
stream nodes. All of the present calculations were performed using 30 
cross-stream nodes. The cross-stream grid spacing varied to include the 
entire width of the flow. 

In order to maintain the program as general as possible, the forward 
step was also varied throughout the flow. The forward step was limited 
so that the quantity of fluid added during the mixing step is a certain 
fraction of the total fluid in the flow to that position in the flow. 
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Table 3 

Definition of Terms in the Generalized GENMIX Program^ 
A = 

® ■ '' 1 * 1 ")'' ^'('e ■ 

C = p u yj./tWj. - 4'j)^o<t>] 

A D 


u a(p^ - p)/p u 

f 0 



^4) = u, f, k, e or g. 
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In the present calculations the forward step was limited so that the 
ratio of additional fluid to the total fluid in the flow was 0.10. This 
ratio was set as low as 0.001 without significant differences in the 
results, however, the larger value reduces the computation time. 

3.5.2 Concentration Fluctuation Model 

The relationships between f and g and y and a in the equations are 
difficult and time consuming to solve. The relationships between these 
quantities were computed once and stored in table form. These tables 
are shown in Appendix D. 

After the solution of the governing equations at e^ch grid node, the 
table was consulted using the specified values of f and g to determine 
values of y and a. 


The probability density function at each grid point is specified by 
the values of y and a. The Dirac delta functions at mixture fractions 
of zero and one can be calculated from Equations (3.8) and (3.9). 


The time averaged temperature, density, and component mass fractions 
can now be calculated at each grid point with the relationship. 


(p = A(pg + + 


a2'ir 


(J)(f) exp 


L- 


£-p 


oy df 


(3.30) 


where <p = T, p, Y^, etc. 

The relationships for 4>(f) provided by the equation of state. The 
integral in Equation (3.30) is evaluated numerically, in a stepwise 
fashion. The value of cp was taken to be constant over several regions 
of the flow. This permitted the error function to be integrated 
independent of the quantity to be averaged. 


CHAPTER IV 

RESULTS AND DISCUSSION OF THE LOCALLY HOMOGENEOUS FLOW MODEL 
4. 1 Introduction 

The major objective of the present study was to systematically evaluate 
a locally homogeneous model of spray evaporation. In order to evaluate 
the model, it was tested against a wide variety of single-phase and two- 
phase flows. This included single-phase constant density jets, single- 
phase variable density jets, two-phase gas-liquid jets and an evaporating 
spray. Theoretical predictions of the axial and radial variations of 
mean and turbulent quantities were compared with the experimental results 
of the present study and results of previous investigations (4, 51, 52, 

53, 54). 

The analytical model used in the present study. Equations (3.12) - (3.19), 
contains various constants which must be specified. Tlie parameters used 
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in the present study were those suggested by Lockwood and Naguib (46) . 

These values are summarized in Table 4. 

This group of constants has been optimized to produce the best overall 
agreement between predictions and measurements for axisymmetric flow. 

One undesirable feature is that and C must be modified in going 
from constant density to variable 2 ^2 density jets (46). This 

is probably due to the rather gross simplifications required to obtain 
model equations for variable density flows, using Reynolds averaging, 
e.g., dropping most density fluctuation terms. The use of Favre averaging 
might eliminate this difficulty, but at the expense of postulating a 
Favre joint probability density distribution. In general, planar jets 
require different constants than axisymmetric jets, which is another 
defect of the turbulence model. Due to this fact, different constants 
are frequently used for the potential core region of axisymmetric jets 
(which is nearly planar, initially), than for the fully developed portion 
of the flow. This modification was not found to be necessary during 
the present study. The same constants were used throughout the flow 
field. 

4 . 2 Rad ial Variation of Mean Quantities 

The theoretical predictions of the radial variation of mean axial velocity, 
mixture fraction, temperature, and liquid mass flux were compared 
with the experimental results found during the present study and several 
previous investigations (4, 51-54). The results were examined in a 
systematic fashion beginning with single-phase constant density jets and 
concluding with the two-phase jets. 

Figure 10 is an illustration of the radial profile of the axial velocity 
for an isothermal single-phase jet, a variable density single-phase jet 
and an evaporating spray. The axial velocity is normalized by the 
velocity at the centerline. The radial coordinate is normalized by 
the axial distance from the injector so that a direct measure of the 
prediction of the width of the flow can be obtained. 

The isothermal single-phase jet is an air jet injected into an air 
environment at the same temperature. The theoretical prediction was 
compared with the experimental results of Wygnanski and Fiedler (51), 
Hetsroni and Sokolov (53), which were obtained with a hot-wire anemometer, 
and the experimental results of the present study, employing the laser 
Doppler anemometer. The experimental results are in the fully developed 
region of the flow extending from 35 injector diameters to 510 injector 
diameters from the injector. The experimental results of the present 
study were in good agreement with the results of the previous air jet 
investigations. The theoretical prediction of the axial velocity profile 
is in good agreement with all of the experimental measurements. 

Analysis of the variable density single-phase jet considers experimental 
measurements within an isothermal sulfur hexafluoride gas jet injected 
into an air environment, obtained during the present investigation. 
Measurements of the velocity were limited to the region far downstream 
of the injector. The theoretical prediction slightly underestimates 
the radial spread of the jet, but is in generally good agreement with 
the experimental results. 
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Table 4 


Constants in the Turbulence Model 


Constant 

Value 

C 

y 

0.09 

c 

1.44 



c = c 

1.89, 1.84^ 



C 

2.8 

®1 



1.0 

a 

1.3 

e 



0.7 

G 

0.7 





a 


Constant density and variable density flows , respectively . 
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Figure 10 Radial Variation of Mean Axial Velocity 
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The effectiveness of the present model for the case of an evaporating 
spray was evaluated by examining a well atomized Freon- 11 spray. The 
experimental results were confined to the region far downstream of the 
injector. Although the model slightly underestimates the spreading rate 
of the spray, the agreement between theory and experiment is relatively 
good. 

The theoretical model was also evaluated on the basis of its ability to 
predict the radial variation of mean mixture fraction. The comparison 
between the predicted and experimental results is shown in Figure 11 
for an isothermal single-phase jet, a variable density single-phase 
jet and an evaporating spray. The mean mixture fraction is normalized 
by the mean mixture fraction at the centerline. The radial coordinate 
is normalized by the axial distance from the injector. 

The ability of the theoretical model to predict the radial variation of 
mean mixture fraction for an isothermal single-phase jet was evaluated 
by comparing the theoretical prediction with experimental results ob- 
tained by Becker, et al., (52). These results were obtained in an 
isothermal air jet.- The theoretical predictions compare well with the 
experimental results. 

The theoretical model also provides a good prediction of the radial 
variation of mean mixture fraction for single-phase variable density 
jets and two-phase evaporating sprays. The experimental measurements for 
both cases were confined to the region far downstream of the injector. 

Both theory and experiment indicate self-preserving flow in this region. 

The locally homogeneous two-phase model was also applied to the two-phase 
air-water jet investigated by Tross (4). Comparison of the radial 
profiles for the axial velocity and mean mixture fraction are shown in 
Figure 12. The model tends to underestimate the radial width of the 
mixture fraction and axial velocity profiles. This disagreement is 
disappointing since the low inertia of bubbles in a continuous liquid 
phase should be favorable to the application of the locally homogeneous 
flow model. In view of the reasonably good predictions for the other 
single-phase and two-phase flows, e.g.. Figures 10 and 11, problems with 
the experiment may be a factor. In particular, the probe used by Tross (4) 
to measure void fraction tends to underestimate the void fraction due 
to surface tension effects. The theory also has unique problems with 
this flow, particularly near the edge o f the jet. As discussed in 
Appendix D, omission of terms involving p’f in the turbulence model is 
particularly questionable in this case due to the very large density 
gradients in the flow. Other terms involving density fluctuations may 
also be more important than in the other flows, involving more modest 
density variations. Finally, the velocity difference between the bubbles 
and the liquid is more appreciable near the edge of the flow, where the 
liquid velocity is relatively low. In view of these factors, the comparison 
between predictions and measurements in Figure 12 is not conclusive, and 
further theoretical and experimental study of this flow would be desirable. 

The radial variation of temperature in an evaporating spray was also 
considered in the present study. Predictions and measurements of the 
normalized temperature decrement for the evaporating Freon-11 spray are 
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Figure 11 


Radial Variation of Mean Mixture Fraction 
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illustrated in Figure 13. Temperature measurements were obtained at 340 
and 510 injector diameters from the injector. Temperature measurements 
could not be conducted closer to the injector, since water vapor in the 
air condensed and froze on the thermocouple junction in this region. 

The radial width of the temperature profile is somewhat overestimated 
by the theory. This is probably because the theory predicts that the 
Freon- 11 should have already evaporated, the lack of any liquid material 
would account for higher temperatures. 

The radial variation of the liquid mass flux is shown in Figure 14. The 
liquid mass flux was measured in the evaporating spray at 3 axial 
positions far from the injector. However, the locally homogeneous theory 
predicts that all of the liquid has vaporized less than 165 inj ector 
diameters from the injector. If the experimental measurements of the 
liquid mass flux are compared with the theoretical prediction of a 
position of 30 injector diameters from the injector, good agreement is 
obtained between the theory and the experimental results. 

Overall, the locally homogeneous theory predicts the radial variation 
of mean quantities reasonably well. In the case of a spray the two- 
phase nature of the flow does not appear to affect the radial variation 
of mean quantities for the present test conditions. Differences between 
the theory and the experiment for an evaporating spray are similar to the 
differences between theory and experiment for a single-phase variable 
density jet. The comparison between theory and experiment for the gas 
liquid jet was the least satisfactory of all the flows studied. However, 
there are unique experimental and theoretical difficulties with this 
flow which deserve further study. 

4.3 Radial Variation of Turbulence Quantities 

The radial variations of the intermittency and the normalized mean 
square mixture fraction for an isothermal air jet are shown in Figures 
15 and 16. The theory is compared to the results of Becker, et al., 

(52) . Fair agreement is achieved for the variation of the intermittency 
prediction over most of the width of the jet. However, the theory 
begins to fill near the edge of the flow. A similar effect is found 
for the variation of the mean square concentration fluctuations. Agree- 
ment is fair near the center of the flow but is less satisfactory near 
the edge of the jet. 

The radial variation of the normalized Reynolds stress is illustrated 
in Figure 17 for an isothermal air jet, a variable density single-phase 
jet, and an evaporating spray. The Reynolds stress is normalized by the 
square of the axial velocity at the centerline. 

The theoretical prediction of the Reynolds stress for an isothermal single- 
phaSe jet is compared with experimental results from Wygnanski and Fiedler 
(51) which exhibit more scatter. In general, the present calculations 
are in fair agreement with the measurements. 
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Figure 13 Radial Variation of Temperature in an Evaporating 
Freon-11 Spray 
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Figure 14 Radial Variation of Liquid Mass Flux in an 
Evaporating Freon-11 Spray 
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Figure 15 Radial Variation of Intermit tency in an Isothermal 
Air-Air Jet [52] 
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Figure 16 Radial Variation of Mean Square Mixture Fraction in 
an Isothermal Air— Air Jet [52] 
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Figure 17 Radial Variation of Reynolds Stress 
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The comparison of the theoretical and experimental results for the 
variable density single-phase jet and the evaporating spray are very 
much alike. Good agreement is obtained between the theory and the exper- 
iment at the larger distances. Deviation between the theory and the 
experiment at the near position is due to experimental error. 

The variation of the turbulent kinetic energy is shown in Figure 18 
for isothermal single-phase jets, an isothermal variable density jet, 
and an evaporating spray. The theoretical model tends to underestimate 
the turbulent kinetic energy for all three flows near the centerline 
of the flow. The experimental data also shows a slight variation in the 
turbulent kinetic energy profiles with x/d which is not predicted by 
the model. 

In general, the theory does a remarkably good job of estimating turbu- 
lence quantities in both the single and two-phase flows. This is 
encouraging, in view of the many simplifications of the turbulence 
equations for variable density flows. Comparing the spray with the other 
flows, it does not appear that the presence of drops had a large influence 
on the turbulence characteristics of the spray. However, it should be 
recalled that the present measurements were limited to the region far 
from the injector, where the spray is very dilute (void fraction greater 
than 99.9%). 

4.4 Axial Variation 

The axial variations of the centerline axial velocity and mixture fraction 
are shown in Figures 19 and 20 for a wide range of single-phase and two- 
phase flows. The axial velocity and mixture fraction at the centerline 
are normalized by their respective values at the injector exit. The 
axial distance is normalized by the injector diameter. In general, 
variation of density ratio between the injected fluid and the surrounding 
gas has a much stronger influence on the axial variation of flow properties 
than the normalized radial variations considered in the previous section. 

The two-phase air-water jet investigated by Tross (4) had the lowest 
density ratio of all the flows examined. In this case, the density ratio 
of the injected fluid to the surrounding fluid was 1 to 1000. The 
predicted decay of the axial velocity and mixture fraction at the center- 
line begins at 1.1 injector diameters from the injector. The theoretical 
prediction of the axial variation of centerline velocity is in good 
agreement with the measurements. However, the theory appears to under- 
estimate the decay of the mixture fraction. The theoretical and 
experimental problems with this flow that were discussed earlier may 
account for these differences. 

The heated air jet analyzed in the present study involved a density 
ratio of the injected fluid to surrounding fluid of 1 to 2. The theory 
predicted that the potential core length would be extended to 5.6 
injector diameters for this flow. The theoretical results concerning 
both the potential core length and the axial decay rates are in good 
agreement with the experimental results obtained by Corrsin, et al., (54). 
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Figure 18 Radial Variation of Turbulent Kinetic Energy 
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The theoretical model was also applied to an isothermal air jet injected 
into air. In this case, the density ratio between the injected and 
ambient fluid is equal to unity. The predicted potential core length 
extends to 7 injector diameters from the injector. The theoretical 
prediction of the velocity decay was in good agreement with the exper- 
imental results of both Wygnanski and Fiedler (51) and the present study. 
Good agreement was also obtained between the theory and experimental 
results of Becker, et al . , (52) for mixture fraction variation. 

A single-phase sulfur-hexafluoride jet that had an initial density 
greater than the surrounding air was also examined. In this case, the 
density ratio of the injected fluid to the surrounding fluid was 
approximately 1 to 0.2. The predictions and measurements of velocity 
and mixture fraction, for this case, are also illustrated in Figures 19 
and 20. The agreement is reasonably good. 

After achieving reasonable agreement between the theoretical predictions 
and the experiments for a wide range of single-phase flows, the locally 
homogeneous model was examined for an evaporating spray. The theory was 
compared with the measurements of velocity and mixture fraction for an 
evaporating Freon- 11 spray conducted in the present study. The results 
are also illustrated in Figures 19 and 20, Agreement between the theory 
and the experiment for the velocity and mixture fraction is significantly 
poorer than for the single-phase jets. In the next chapter this difference 
will be shown to be largely the effect of slip and loss of thermodynamic 
equilibrium. 

The effect of the loss of kinematic and thermodynamic equilibrium can 
also be observed in the axial variation of temperature in an evaporating 
spray. Figure 21 is an illustration of the comparison between the pre- 
diction of the axial variation of the normalized temperature decrement 
with the experimental results. The temperature decrement (the difference 
between the centerline and ambient temperatures) is normalized by its 
value at the injector exit. The theoretical prediction increases to a 
maximum value of approximately 1.5 and then decreases throughout the rest 
of the flow. The increase, in the temperature decrement is predicted by 
the equation of state. The maximum value of the temperature decrement 
corresponds to the position where the locally homogeneous model predicts 
that all of the liquid has vaporized. In contrast, the measurements 
suggest a much slower development of the mean temperature, and liquid was 
present at all positions in the flow. 

Table 5 lists the variation of liquid mass flux along the jet centerline 
for the evaporating Freon- 11 spray. The theoretical model predicted that 
all of the liquid has vaporized beyond 30 injector diameters from the 
inj ector . 

The ability of the model to predict the axial variation of the concentra- 
tion fluctuations was also examined. Figure 22 is an illustration of the 
theoretical prediction of g and the experimental results of Becker, et al., 
(52). Reasonable agreement was obtained between the theory and the 
experiment, at least for this cc5hstant density single phase flow. 
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Table 5 


Axial Variation of Liquid Mass Flux for an 
Evaporating Freonr-11 Spray 
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.011 






Figure 22 Axial Variation of Concentration Fluctuations in an Isothermal Air-Air Jet 
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4 . 5 Droplet Size 

The axial and radial variation of the droplet size is shown in Table 
6. The Sauter mean diameter remains nearly constant across a radial 
profile and increases slightly with axial distance. The increase in 
Sauter mean diameter with distance is caused by the rapid evaporation 
of small drops leaving only the larger drops at the downstream position. 
Naturally, features of this type cannot be considered using a locally 
homogeneous flow model . 


CHAPTER V 

DROPLET LIFE-HISTORY CALCULATIONS 
5. 1 Theoretical Model of Droplet Life-History 

The lifetime of a single drop evaporating in a convective environment 
was modeled in the present study. The calculations considered a drop 
in the evaporating spray considered in the present test program. Using 
these results, a direct assessment of the magnitude of the slip and loss 
of thermodynamic equilibrium between the phases could be obtained. In 
completing the calculations, an estimation of the local environment of a 
drop within the spray was obtained from the results of the locally 
homogeneous flow model . 

The drop life-history calculations were calibrated against measurements of 
a single supported drop evaporating in an air stream of known velocity 
and temperature, at atmospheric pressure. In this manner, errors in the 
drop life-history calculations (due to uncertainties in selecting proper 
average properties for transport parameters) were minimized. 

Throughout the analysis, a single droplet in a dilute spray is considered. 
The droplet is assumed to be spherical in shape and unaffected by the 
presence of any other droplets. The gas phase is taken to be a quasi- 
steady continuum, fully adjusted to the steady state structure for the 
imposed droplet boundary conditions at each instant of time. The pressure 
throughout the flow is assumed to be constant. 

The behavior of all species in the gas phase is approximated by the ideal 
gas law. The surrounding gas is assumed to have negligible solubility 
in the gas phase. The partial pressure of the vaporized liquid is given 
by Equation (B.19). 

The effects of chemical reaction are assumed to be negligible and the 
Lewis number is taken to be unity. The other fluid properties are 
evaluated at the local conditions in the flow. The .values of the fluid 
properties and a description of how the properties are evaluated are 
presented in Appendix E. 

Under these assumptions the governing equations are: 
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Table 6 


Variation of Sauter Mean Diameter (]Jm) in 



Evaporating 

Freon-11 

Spray 


X 


r/x 


d 

.026 

.050 

.075 

170 

27.2 

27.6 

27.6 


340 


29 .4 


29.5 


29.2 
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Conservation of Mass 
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dx 
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Conservation of Momentum 
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(5.2) 


Conservation of Energy 
dT 12ANj^ 

ji* = 5 ln(l + B ) [C (T - T ) - B h. ] (5.3) 
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^f Re 6 
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(5.4) 


Re = 


D(Up - u) 


(5.5) 


1/2 1/3 4/3 

= 1 + 0.276 Re ' Pr ^ /(I + 1 . 232/ (RePr) ^ ) (5.6) 


K = (Xn - Y„ )/(l - Y„ ) 

y Fgp Foo^ Fgp^ 


(5.7) 


The drop position is obtained from the relationship 
t 


X = u dt 

P 

0 


(5.8) 


The drops and the gas phase are assumed to be at the same velocity and 
temperature throughout the potential core. 


For the spray, the ambient conditions of the drop were obtained as a 
function of position from the locally homogeneous flow model. For the 
calibration experiment, the ambient state was defined by the predictions 
of the locally homogeneous model. In this case. Equations (5.1) and 
(5.3) were transformed to a function of time by noting 


d ^ d 

^p dx dt 


(5.9) 
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Equations (5.1) - (5.8), along with the vapor pressure relationship 
between Yp p and T , obtained from Equation (B.17), are sufficient to 
solve for the temperature, diameter, position and velocity of the 
droplet. The system of equations was integrated on a digital computer, 
using a fourth-order predictor-corrector method. 

5 . 2 Calibration Apparatus 

The purpose of the experimental apparatus was to provide a convective 
environment for the observation of an evaporating drop. These exper- 
imental results were used to calibrate the drop life-history calculations. 

A schematic of the overall experimental apparatus is shown in Figure 
23. Motion pictures of the droplet evaporation process were taken to 
record the variation of droplet size. The droplet was suspended on a 
fine wire thermocouple which also served to measure the droplet temper- 
ature. A Pitot-static velocity probe was used to measure the velocity 
of the surrounding gas. 

Turbulent air flow at the exit of a pipe provided the convective flow 
field for the droplet evaporation process. Fully developed turbulent 
pipe flow was obtained from a pipe 32 mm diameter and 610 mm in length. 
Four 16 mesh screens were placed near the pipe inlet to assist in deve- 
loping the flow. 

The air flow rate was metered and controlled with a combination of a 
pressure regulator and critical-flow orifice. The dry compressed air was 
supplied by Matheson Gas Products. A Matheson model 4 regulator was 
used to control the air pressure upstream of the critical flow orifice. 

The upstream pressure was measured with a Heise model C-54848 pressure 
gauge. The orifice was calibrated with a Precision Scientific wet-test 
meter (.283 5,/rev.). 

The drop temperature was determined by suspending the drop on the 
junction of a fine wire thermocouple. The thermocouple was manufactured 
from 0.076 mm chromel-alumel wire. The junction had a diameter 
of 0.5 mm. The reference junction for the thermocouple was placed in an 
ice bath. A Hewlett-Packard X-Y plotter, model 7044A, was used to record 
the droplet temperature variation with time. 

Motion pictures of the drop evaporation process were used to measure 
the variation of droplet size with time. A Photosonic 1-B motion 
picture camera driven with a Kepco SM 36-5 AM d.c. power supply were 
used throughout the investigation. The film speed was indicated by a 
timing light on the camera activated by an Adtrol Electronics pulse 
generator, model 501. An event marker was used to correlate the photo- 
graphic and temperature measurements. Kodak plus-x reversal film was 
used for all tests. 

Backlighting for the photographic measurements was supplied by a Pek, 
model 401A, arc lamp using a 75 watt mercury bulb. The light was focused 
into a parallel beam using the optics located in the arc lamp. A diffuser 
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Figure 23 Schematic Diagram of the Drop-Life History Calibration Apparatus 
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screen was employed to equalize the background intensity of the light 
beam. 

The air velocity was determined with a Pitot-static probe. The probe 
was connected to a Meriam model 34FB2 micromanometer in order to measure 
the difference between the static and stagnation pressure. 

Measurements of the droplet evaporation process were controlled by a 
mechanical timer. As soon as a stable droplet was established on the 
thermocouple bead, the timer was activated. At the beginning of the 
cycle, the electrically driven motion picture camera was started and 
allowed to reach operating speed. The event marker and X-Y plotter 
were then activated simultaneously. The process was allowed to continue 
for approximately 4 seconds before the camera and the X-Y plotter were 
turned off to complete the test. 

The temperature trace consisted of three regions; first, a region of 
rapidly decreasing temperature; second, a constant temperature region; and 
third, a region where the temperature increases back to the ambient gas 
temperature. The initial temperature decrease represents the droplet 
approaching an equilibrium condition with the moving air stream. The 
constant temperature is the wet-bulb temperature for the evaporation 
process. At this condition, all of the energy transferred to the drop 
is used to vaporize the liquid. The final region follows the total 
evaporation at the liquid when the thermocouple returns to the ambient 
gas temperature. 

The analysis of the droplet evaporation process was confined to the constant 
temperature region. The films were analyzed on a frame by frame basis 
using a Vanguard motion picture analyzer. The motion picture analyzer was 
used to view film and measure the size of the droplet in each frame. 

5.3 Calibration Test Results 

The purpose of the calibration tests was to establish the proper 
reference condition to evaluate the fluid properties in Equations (5.1) 
to (5.7). The calibration experiment provided a well-defined test 
condition to evaluate the drop evaporation model. The ambient velocity, 
concentration and temperature were constant with time. The theoretical 
and experimental results for the variation of droplet size and temperature 
were examined. 

Several test conditions were considered with ambient velocities from 
3 m/sec to 6 m/sec. This velocity range was selected so that the drop 
Reynolds Numbers in the calibration experiment spanned the range encountered 
in the spray experiment. The ambient gas temperatures were in the range 
18-22 C. 

The closest agreement between predictions and the measurements was 
achieved when gas phase properties in Equations (5.1) - (5.8) were 
evaluated at the droplet surface temperature. Figure 24 is an illustration 
of some typical theoretical and experimental results. The drop diameter 
variation is plotted as a function of time for two different gas velocities. 
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Figure 24 Drop Life-History Calibration Results 
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The agreement between prediction and measurements is seen to be quite 
good. The wet-bulb temperature of the drop is not a strong function of 
drop size or velocity. The predicted value is within 3 C of the measured 
value, which is satisfactory for present purposes. 

Based on the results of the calibration experiments, the drop surface 
temperature was employed as the reference temperature for properties 
in the present drop life-history calculations for spray conditions. No 
universality of this reference condition should be implied. The uncer- 
tainties in average properties and the transport correlations are such 
that calibration should be undertaken for each system to be investigated (33) 

5.4 Droplet Life-History Calculations in Evaporating Spray 


The comparison between theory and experiment for the spray suggests that 
slip and loss of thermodynamic equilibrium are retarding the development 
of the flow. This is not completely convincing, however, since a large 
number of assigned constants could be adjusted to provide better agree- 
ment between theory and experiment. In particular, the values of and 
C for variable density flow are somewhat arbitrary (46) and the 2 
^2 entrainment rate of the flow is quite sensitive to these parameters. 

In order to provide a more direct evaluation of the locally homogeneous 
assumption, drop life-histories were computed in the actual spray 
evaporation process. This involves determining the variation in drop 
size, velocity, and temperature as a function of position in the flow. 

The drop life-history calculations conducted in the present study considered 
a drop passing along the centerline of the spray. As an estimate of the 
environment of the drop, gas velocities, temperatures and composition 
were obtained from the locally homogeneous flow prediction. The drop 
is assumed to be in thermal and kinematic equilibrium with the gas 
phase, until it reaches the end of the potential core. 

The calculations were conducted for the conditions of the evaporating 
spray. Results for drops having initial diameters of 10, 30 and 50 pm 
are illustrated in Figure 25. The diameter, velocity and temperature of 
the drop, and the velocity and temperature of the gas at the centerline 
of the spray are given as a function of the distance from the injector. 

The drop velocities decay much more slowly than the gas velocities. 
Appreciable levels of slip occur in the region just downstream of the 
potential core. The difference between the drop velocity and the gas 
velocity indicates the magnitude of slip present. For these calculations, 
where the rate of development of the flow is overestimated, the gas 
temperature falls rapidly near the injector. Farther downstream the 
gas temperature rises toward the ambient temperature, while the drop 
temperatures tend toward a constant value, which is below the gas tempera- 
ture. 

The drop life-history calculations indicate that drops are present for 
x/d^ greater than 200, which agrees with experimental observations. These 
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Figure 25 Drop Life-Histories in an Evaporating Freon— 11 
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estimations are still conservative, since gas properties were taken 
from the locally homogeneous flow model, which overestimates the rate 
at which jet conditions approach the ambiance. 

The results shown in Figure 25 indicate that slip and loss of thermo- 
dynamic equilibrium are important for the present test conditions. This 
accounts for the overestimation of the rate of development of the spray 
using the locally homogeneous flow model, even for drops as small as 
10 pm. From this evidence and the results of Khalil and IVhitelaw (18) 
and Shearer and Faeth (20), it appears that slip must be considered unless 
information is only required far from the injector. 


CHAPTER VI 

SUMMARY AND CONCLUSIONS 


6. 1 Summary 

The overall objective of the present study was to develop and evaluate a 
locally homogeneous flow model of spray evaporation which employs a 
second-order turbulence model. The locally homogeneous flow assumption 
implies no velocity difference and thermodynamic equilibrium between 
the phases. The turbulence was represented by a k-e-g turbulence model 
employing a clipped Gaussian probability density function for mixture 
fraction fluctuations. 

The theoretical model of the spray process was based upon the Reynolds- 
averaged form of the conservation equations. A modified version of 
GENMIX (45) combined with an appropriate equation of state was used to 
solve the governing equations. The program could be applied to both 
single-phase and two-phase flows. 

Experimental measurements of mean and fluctuating velocities were 
conducted in an isothermal air-air jet, an isothermal sulfur hexafluoride- 
air jet, and an Evaporating Freon- 11 spray. A single-channel LDA was 
used to conduct these measurements and the results were combined to 
determine the Reynolds stress and turbulent kinetic energy. The measure- 
ments were conducted from 170 to 510 injector diameters from the injector. 

Experimental measurements of the mixture fraction were conducted in 
the sulfur hexafluoride-air jet and the evaporating spray. The flows 
were isokinetically sampled, based on the velocity measurement at each 
location. The mixture fraction measurements were extended to the region 
far downstream of the injector (170 - 510 injector diameters). 

The Freon- 11 spray was further investigated by measuring the gas tem- 
perature and drop size throughout the spray. The temperature was 
measured at positions from 200 to 510 injector diameters from the injector. 
Water vapor condensing and freezing on the thermocouple junction pro- 
hibited measurements at closer positions to the injector. The drop size 
was measured using a magnesium oxide impaction technique at positions of 
170 to 510 injector diameters from the injector. 
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The experimental results of the present study, combined with the 
experimental results of previous studies (4, 51-54), consider a wide 
range of flow conditions. The theoretical model was systematically 
evaluated against each of these flow conditions. 

In order to evaluate the possible effect of slip in the two-phase 
evaporating spray, the behavior of a single drop in the spray was 
examined. A single drop was assumed to follow the centerline of the 
spray and the ambient conditions of the drop were those predicted by 
the locally homogeneous flow model. The axial variation of temperature, 
velocity and size provided an indication of the degree of slip and the 
loss of thermodynamic equilibrium between the two phases. 

6. 2 Conclusions 

The major conclusions and observations of this study are: 

1. The locally homogeneous flow assumption in conjunction with a 
second-order turbulence model provides a good prediction of the 
axial and radial variation of mean quantities in single-phase 

jets. Fair agreement between the experimental and analytical results 
was achieved for the variation of turbulence quantities. The 
present model tends to underestimate the radial width of some 
turbulence quantities. 

2. The locally homogeneous flow model provides a reasonably good 
estimate of the radial variations of mean quantities in an evapora- 
ting spray. However, the locally homogeneous flow model tends to 
overestimate the rate of development of the flow, even for a spray 
having a Sauter mean diameter of 30 ym. This is illustrated by 
examining axial variation of velocity and mixture fraction. While 
the rate of reduction of velocity and mixture fraction defect is 
predicted correctly, the measured positions where a particular value 
is reached lag the predictions by 15 and 40 percent, respectively. 

3. The locally homogeneous flow model can be an extremely useful design 
tool. Because only a minimum amount of initial information is 
required, the locally homogeneous model can be easily applied to a 
wide range of flows to provide a reasonable initial estimate of flow 
characteristics . 

4. Analysis of the lifetime of a single drop shows that slip and loss 

of thermodynamic equilibrium are important even for the well -atomized 
spray of the present study. It appears that a Sauter mean diameter 
less than 10 ym would be required for quantitative accuracy with the 
locally homogeneous flow model. These small sizes are needed due to 
relatively fast decay of the flow due to the small injector diameter. 
Since injector diameters of sprays are generally small, this finding 
implies that accurate analysis of most practical sprays will require 
the use of models which allows for slip and finite rate processes 
between the phases. 
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6. 3 Recommendations for Further Study 

Although the locally homogeneous flow model tends to overestimate the 
rate of development of the flow, it provides a useful qualitative picture 
of the spray evaporation process. The locally homogeneous model should 
be extended to consider two-phase combusting sprays. It is expected 
that slip and loss of thermodynamic equilibrium would be important in 
combusting flows, however, the simplicity of the approach and the lack 
of any need for detailed injection characteristics would establish the 
model as a useful design tool. 

A more complete examination of the spray evaporation and combustion process 
will require the use of a two-phase flow model. Systematic validation of 
current two -phase flow models will be required. Further work must 
determine the effect of particles on turbulence production dissipation, 
and scale, the effect of particle diffusion in sprays, and the limit 
of the dilute spray approximation. 

A standard procedure for defining injector characteristics must also 
be developed. More complex slip models will require very detailed 
information specifying the initial conditions (initial drop size 
distribution, initial velocity distribution, etc.). Unless this infor- 
mation can be accurately determined, it is doubtful that the two-phase 
slip models will be any more accurate than a locally homogeneous model. 
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APPENDIX A 

EVALUATION OF TURBULENCE QUANTITIES 

Consider a single particle moving through the fringe pattern shown on Fig- 
ure 26(A) produced by the beams crossing. The distance between fringes is 

d^ = X/(2 sin C(J)/2)) (A.l) 


The velocity is 

u = Fd^ (A. 2) 

The measured frequency is directly proportional to the output of the 
single processor. 

F = KE (A. 3) 


Combining Equations (A.l), (A. 2), and (A. 3) yields 


E = Cu 


(A. 4) 


where 


C = 2 sin (cl)/^/2)/KX (A. 5) 

In turbulent flows the velocity consists of mean and fluctuating com- 
ponents producing mean and fluctuating components of the output voltage 

(E + e’) = (u + u')C (A. 6) 


However, it is not required that the flow be perpendicular to the fringe 
pattern. This is illustrated in Figure 26(B). 

E^ + e' = C[u sin a + v cos a + u' sin a + v' cos a] (A. 7) 


The average of the output voltage is determined by averaging Equation (A. 7) 


E = C[u sin a + V cos a] (A* 8) 

The relationship between the fluctuations output voltage and the tur- 
bulence quantities is obtained by squaring Equations (A. 7) and (A. 8) and 
combining the resulting equations 




Figure 26 Sketch of a Particle in the LDA Fringe Pattern 
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,2 ,2 . 2 

' = C lu’ sin a + 


2 2 


C [u’ sin a + V* cos a + 2u'v’ cos a sin a] 


I A. 9) 


2 2 

Values for u' , v’ , and u*v* may be derived by 3 measurements at 3 
different values of a. 


For the special cases = 45® and = 135®, u'v’ may be determined 
from the following equari“ 


:ion 


u'v' = - e'p/2C 


(A. 10) 


APPENDIX B 
EQUATIONS OF STATE 


B. 1 Isothermal Air Jet 

Experimental results from the present study as well as earlier measure- 
ments [51-53] were considered in these calculations. The test conditions 
involved air flowing into stagnant air at room temperature and pressure 
(298 K, 101 kPa). 

Since the flow has a constant temperature, pressure and composition 

T = T CB.l) 


P = P» = Po,/Ra'^» (B-2) 


B, 2 Heated Air Jet in Air 

Experimental results reported by Corrsin [54] were considered for this 
flow. The test conditions involved air heated to 300 C, flowing into 
stagnant room air (298 K, 101 kPa) . 

The air was assumed to have a constant specific heat. Therefore, from 
Equation (3.4), the enthalpy of the gas mixture is 

= Cp + (1 - f) (T„ - T^)] (B.3) 

A 

where the reference temperature, Tj^, was 298 K. A second expression for 
the mixture enthalpy can be obtained from Equation (3.3), as follows 


h = Cp^CT - 


(B.4) 
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Given f, and the initial and boundary conditions. Equations (B.3) and 
(B.4) provide the mixture temperature. 

The density o£ the mixture was obtained from the ideal gas law 

P = P/Ry^T (B .55 


B i.3 Dense Gas Jet in Air 

The experimental results for the dense gas jet in air were obtained during 
the present investigation. The flow was isothermal, involving pure sulfur 
hexafluoride injected into stagnant room air (298 K, 101 kPa) . Noting 
that the mass fraction of sulfur hexafluoride is unity at the injector 
exit and zero in the ambiance. Equation (3.1) yields 

Yg = f (B.6) 

= 1 - £ (B.75 


The density o£ the mixture can be obtained from Equation (3.5) 


P = [f/Pg + (1 - f)/p^] 


(B.8) 


where the density of the individual components is calculated from the 
ideal gas law at the temperature and total pressure of the mixture. 

B.4 Air Jet in Water 

Data obtained by Tross [4] was compared with the analysis for this flow. 
The experiments involved a pure air jet injected into a stagnant bath of 
water at 298 K, 101 kPa. While expansion of the gas through the injector 
resulted in a gas temperature below room temperature at the injector exit, 
this enthalpy defect is negligible due to the large heat capacity of the 
water. Therefore, the flow was assumed to be isothermal. The vapor 
pressure of water in the air was also neglected. 

Under these assumptions, the mass fractions of air and water can be found 
from Equation (3.1) as follows 

= f (B.9) 

= 1 - f (B.IO) 


The density of the flow is given by Equation (3.5) 
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p = [f/p^ + (1 - f)p^] 


(B.ll) 


The density of both the water and the air was assumed to be constant . 

The air density was computed assuming an ideal gas . 

B.5 Evaporating Spray in Air 

The data for an evaporating spray was obtained in the present investi- 
gation. An air atomizing injector was used with Freon-11 as the evapor- 
ating material. Therefore, both air and Freon-11 leave the injector. 

The spray was injected into stagnant room air at 298 K, 101 kPa. 

Given the mixture fraction, the mass fractions of Freon-11 and air can 
be determined from Equation (3.1) 


Y 


F 


^FO^ 


(B.12) 


Y 


A 



(B.13) 


where Y is known from the injector operating condition, 
r U 

The enthalpy of the mixture can be determined from Equation (3.3) 


h = h^ f + h. (1 - f) 
0 A<»' ^ 


(B.14) 


The enthalpy of air is given by Equation (B.4) . The enthalpy of Freon-11 
in the liquid and gas phases were expressed as 


•^Ff = Sp/" - V 


h = C (T - T ) h 

Ppg *^fgR 


(B.15) 

(B.16) 


where T^ was taken to be 298K. 

K 

The enthalpy at the injector exit is the sum of the contributions due to 
Freon- 11 and the injector air. 


h 


0 


^Ao\o ^FO^FO 


(B.17) 


The gas phase consists of both Freon-11 vapor and air. Assuming ideal gas 
behavior, the total pressure is equal to the sum of the partial pressures 
of the individual components. 
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= Ep. = 


Fg 


(B.18) 


When the liquid is present, the partial pressure of Freon-11 in the gas 
phase must equal the vapor pressure. The vapor pressure of Freon- 11 
can be correlated by the following expression 

logio Ppg = A' - BVT (B.19) 


Then, since mole fractions in the gas phase are proportional to the 
partial pressures, the mass fraction of Freon-11 in the gas phase is 

V ^ Vpg VPa^ 

Then 


Y = Y - Y 
Ff F Fg 


The enthalpy of the mixture can also be determined from Equation 


^ ^A^Fg ^Ff^Ff ^Fg^Fg 


(B.21) 

(3.4) 

(B.22) 


As long as the mixture is saturated, the mixture temperature and com- 
position can be calculated at any mixture fraction using Equations 
(B. 12) - (B . 22) . A temperature is estimated and the partial pressure of 
Freon vapor is determined from Equation (B.19). The mass fractions and 
mixture enthalpy are then calculated from Equations (B . 20) - (B . 22) . The 
enthalpy computed in this manner is compared with the enthalpy determined 
directly from Equation (B.14). An interval halving procedure was em- 
ployed to adjust the temperature until both calculated values of the 
mixture enthalpy agreed. 


Once all the liquid has vaporized, 
and Freon-11 vapor can be determine 


Y ^ = 0 and the mass fractions of air 
d^airectly 



(B.23) 


Y 


A 



(B.24) 


The mixture temperature can then be determined from Equations (B.4), 
(B.14) , (B.16) and (B.22) . 

The mixture density is determined from 



The density of the Freon-11 liquid was taken to be a linear function of 
temperature. The air and vapor densities were obtained from the ideal 
gas expression at the temperature and total pressure of the mixture. 

B. 6 Property Data 


The property data [63] employed in the calculations described in 
Sections B.1-B.5 are summarized in Table 7. 
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Table 7 


Summary of Properties Used in the Equation 
of State Calculations® 


Air: 

M = 28.966 kg/kg mole 

C = 1.005 kJ/kg K 
P 

Sulfur Hexafluoride: 

M = 146.05 kg/kg mole 
Water : 

p = 997.1 kg/m 
Freon-11 : 

M = 137.37 kg/kg mole 

C = 0.674 kJ/kg K 
Pg ^ 

C , = 0.879 kJ/kg K 
Pf 

= 171.17 kJ/kg 
= 2143.7 - 2.235T kg/m^ 
log^QPg = A’ - B’/T kPa 
A’ = 6,7828 
B’ = 1416.1 

in K, Universal gas constant = 8.3143 kJ/kg mole K, = 298 K, 
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APPENDIX C 

PROBABILITY DENSITY FUNCTION TABLE 

Evaluation of y and a for P(f) 

The values of o and y are given as a function of f and g in 
Tables 8 and 9. 
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Table 8 


C 7 as a Function of f and f/g 


1/2 



f =.001 

.01 

.05 

.1 

.2 

.3 

,4 

.5 

1.0 

3.0 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

2.5 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

1 . 

2.0 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.02 

1.04 

1 . 

1.5 

1.09 

1.09 

1.09 

1.09 

1.09 

1.09 

1.12 

1.27 

1 . 

1.2 

1.17 

1.17 

1.17 

1.17 

1.17 

1,17 

1.34 

1.92 

1 . 

1.1 

1.22 

1.22 

1.22 

1.22 

1.22 

1.25 

1.55 

2.90 

10 , 

1.0 

1.29 

1.29 

1.29 

1.29 

1.29 

1.37 

1.99 

10 . 

20 . 

0.9 

1.38 

1.38 

1.38 

1.38 

1.38 

1.59 

3.50 

20 . 

30 . 

.85 

1.43 

1.43 

1.43 

1.43 

1.43 

1.80 

6.40 

30 . 

40 . 

.80 

1.50 

1.50 

1.50 

1.50 

1.53 

2.15 

10 . 

40 . 

50 . 

.70 

1.67 

1.67 

1.67 

1.67 

1.87 

4.70 

20 . 

50 . 

60 . 

.60 

1.90 

1.90 

1.90 

1.93 

2.85 

10 . 

30 . 

60 . 

70 . 

.52 

2.17 

2.17 

2.17 

2.27 

7.80 

20 . 

40 . 

70 . 

80 . 

.40 

2.90 

2.90 

2.90 

4.0 

10 . 

30 . 

50 . 

80 . 

90 . 

.35 

3.40 

3.40 

3.50 

12.0 

20 . 

40 . 

60 . 

90 . 

100 . 

.30 

4.10 

4.10 

4.60 

20 . 

30 . 

50 . 

70 . 

100 . 

110 . 

.25 

5.10 

5.10 

10.7 

30 . 

50 . 

60 . 

80 . 

110 . 

120 . 

.20 

6.65 

6.65 

20 . 

40 . 

50 . 

70 . 

90 . 

120 . 

130 . 

.15 

9.40 

9.90 

30 . 

50 . 

60 . 

80 . 

100 . 

130 . 

140 . 

.12 

12.3 

17.1 

40 . 

60 . 

70 . 

90 . 

110 . 

140 . 

150 . 

.11 

13.6 

28.5 

50 . 

70 . 

80 . 

100 . 

120 . 

150 . 

160 . 

.10 

15.2 

40 . 

60 . 

80 . 

90 . 

110 . 

130 . 

160 . 

170 . 
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Table 9 


y as a Function of f and Cf/f 



a/f 




£ 



.1 

.3 

.4 

.45 

.48 

.5 1.0 


.4 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

.5 

.996 

.996 

1. 

1. 

1. 

1. 

1. 

.8 

,955 

.995 

.980 

.975 

.985 

1. 

1. 

1.0 

.901 

.905 

.935 

.960 

.975 

1. 

1. 

1.25 

.810 

.830 

.885 

.930 

.965 

1. 

1. 

1.5 

.675 

.730 

.830 

.905 

.955 

1. 

1. 

2.0 

.376 

.490 

.710 

.845 

.930 

1. 

1. 

2.5 

.005 

.255 

.590 

.785 

.905 

1. 

1. 

3.0 

- .430 

0.00 

.450 

.720 

.885 

1. 

1. 

3.5 

- .92 

- .26 

.305 

.665 

.860 

1. 

1. 

4.0 

- 1.38 

- .55 

.180 

.600 

.835 

1. 

1. 

5.0 

- 2.46 

- 1.10 

- .03 

.48 

.785 

1. 

1. 

6.25 

- 3.92 

- 1.70 

- .36 

.335 

.725 

1. 

1. 

7.14 

- 5.05 

- 2.18 

- .59 

.23 

.68 

1. 

1. 

8.33 

- 6.64 

- 2.83 

- ,86 

.09 

.62 

1. 

1. 

8.99 

- 7.12 

- 3.10 

- 1.0 

.0 

.59 

1. 

1. 

10.0 

-9.02 

- 3.60 

- 1.29 

- 1.28 

- .14 

1. 

1. 

20.0 

-25.1 

-10. 

-10. 

-10. 

-10. 

1. 

1. 

37.5 

-46.7 

-10. 

-10. 

-10. 

-10. 

1. 

1. 
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APPENDIX D 

ORDER OF MAGNITUDE ANALYSIS OF GOVERNING EQUATIONS 

The order of magnitude analysis considered the equations of conservation 
of mass, momentum, mixture fraction, turbulent kinetic energy, turbulent 
dissipation, and concentration fluctuations. The governing equations for 
this analysis follow those proposed by Gosman, Lockwood, and Syed [47] . 
These equations are limited to steady, axisymmetric, compressible flow. 
The order of magnitude analysis will not be conducted for dimensionless 
equations in order to preserve the notation used in the calculations. 

For boundary layer flows, where the thickness is very small compared to 
the axial distance, it is assumed 

u,x,p ^1; r 'v. 6 (D.l) 

The Reynolds averaged equation of conservation of mass for a steady 
compressible flow is 

(a) (b) (c) (d) 

It * It * ? IF * F IF = ° 


The Boussinesq approximation is employed to model the density-velocity 
correlations appearing in terms (b) and (d) in equation (D.2). 


P'u' 


"^t 9p 

a dx 
P 


(D.3) 


p'v' 


o 8r 
P 


(D.4) 


Substituting Equations (0.3) and (D.4) into Equation (D.2), and applying 
the orders of magnitude indicated in Equation (D.l) yields 


(a) 


(b) 


(c) 


(d) 


dx 


■(pu) 


9x 


(Ut- 


3p 


9x 


-) + 


8r 


(rpv) - 


8r 


•Cw^- 


3p 


9r 


-) = 0 (D, 




V /6 




1 
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where the order of magnitude has been indicated below each term. In 
order that Equation (D.5) not be trivial, terms (c) and (d) must be 
retained. This implies that 

V 5 , 5^ (D.6) 

With this assessment, term (b) in Equation (D.5) is higher-order than 
the rest. Therefore, Equation (D.2) may be rewritten as 

-|r ^ Crpv“) = 0 CD. 7) 


where 


pv® = pv + p*v' 


CD. 8) 


The combination of mean and fluctuating terms given by Equation (D.8) 
provides a convenient formulation for the remainder of the analysis. 
Razdan (64) applies a similar substitution for a two-dimensional planar 
boundary layer. 


The momentum equation in the axial direction may be expressed 


(a) (b) Cc) (d) (e) 

pu + pv = a(p^ - P) - CPu'u') - (rpu’v') 

(f) (g) (h) Ci) 

^ (up'u') i Crvp'v') ^CP'u’u') ^ |^(rp?u’v’) 


(j) (k) 


- p'u' 


9u 

9x 



(D.9) 


when all terms involving the laminar viscosity have been neglected, which 
is well- justified for a turbulent free-jet flow (47). 


As before, the Boussinesq approximation. Equations (D. 3) and C^-^^), can 
be used to model terms (f), (g) , (j) and (k) . Gossman, et al., (47) 
suggest the following expressions for the Reynolds stress terms: 


pu'u' = -2y^ 


3u ^ 2 ^t - 9p 

9x 3 ^ ^ 9x 


(D.IO) 


pu' v' 


f 9u 



+ 


(D.ll) 
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The order of magnitude of pu'u' is indicated by the first term on the 
right hand side of Equation (D.IO). Therefore, only this term will be 
considered in the following in order to simplify the formulation. From 
Equations (D.l) and (D.6), the second term on the right hand side of 
Equation (D.ll) is higher-order than the first and can be dropped 
immediately. The triple correlation terms (h) and (i) in Equation (D.9) 
present special problems. These terms may not be small. However, since 
suitable correlations have not yet been developed for them, they are 
neglected, similar to Gosman, et al . , (47). 

From this discussion. Equation (D.9) may be written, with the order of 
magnitude of the terms indicated, as follows: 


pu 


(a) 
3u 


(b) 


-- 9u 


9x 


+ pv 


9r 


(c) (d) 

= a(p^ - 3) + 2 


9 , 9u 

- ^^t 


(e) 


(rp, 


(f) 


9u 


r 9r ^ t 9r 
1 

(j) 


-) + 


9x 


(u 


(k) 


^t 9p 9u ^t 9p 9u 
+ ^ + 


9x 9x 

.2 


a 9r 9r 
P 


(g) 

iL.) + J_ -3_ 

9x ^ r 9r '' a 9r 


(D.12) 


Due to their higher order, terms (d) , (f ) , (g) and (j) may be eliminated. 
With these simplifications. Equation (D.9) becomes 


pu 


9u 

9x 


+ pv° 


9u 


1 ^ 

r 9r 


(ry. 



a(P«, - P) 


(D.13) 


The conservation of mixture fraction equation can be written as follows: 
(a) (b) (c) (d) (e) 


pu 


9f 

9x 


+ pv 


9f 

9r 


9x 


(f) 


9r 


(rvp'f’) - 


9x 


- (Pu'f) - 

Cg) 

(p’u’f ') - - 


-g— (rpv'f) ^(up’f) 

(h) (i) 


9_ 

'9r 


(rp’v’f’) - p'u' 


9f 

9x 


9f 


(j) 


(D.14) 
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Similar to conservation of momentum, the triple correlation terms 
(g) and (h) will be neglected (47). 


Terms (e) and (f) contain the scalar correlation p'f*. This correlation 
could be modeled statistically. Since p is solely a function of f 
under the assumptions of the present analysis, the time averaged value 
of the correlation can be determined from the relationship 

1 

I p,f, df (D.15) 

0 


Equation (D.15) can be written, using the basic Reynolds averaging 
definitions, as follows 


p»f ' 


1 

(p - p)(f - f) P(f) df 

0 


(D.16) 


Under the assumptions of the present analysis, the equation of state 
provides p as a function of f. The solution of the flow equations 
provides p, T and a and p (which are required to complete the specifica- 
tion of P(f). With this i nform ation the integral in Equation (D.16) 
can be evaluated to yield p’f, and thus terms (e) and (f) in Equation 
(D.14). 

If the magnitude of the fluctuations is not large. Equation (D.15) 
can be approximated. Neglecting higher order terms, since p = p(f), 
we have 


P - P = Cf - f) CD. 17) 

where the derivative is to be evaluated at f. Substituting Equation 
(D.17) into (D.16) yields 

1 

= J Cf - P(f)df CD. 18) 

0 


Equations (D.16) and (D.18) reveal that p’f depends very strongly on 
the form of the probability density function used in the model. Since 
the probability density function is not kno wn ve ry accurately, Gosman, 
Lockwood and Syed (47) chose to ignore the p’f terms. This model gave 
reasonably good predictions of experimental results; therefore, the same 
procedure was chosen in the present study. This introduces no additional 
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uncertainty for most of the flows considered in this study, since the 
density vairations are no greater than those encountered in Reference 
(47). The exception is the air jet into water, where the equation of 
state illustrated in Figure 7 indicates very large values of dp/df in the 
regio n near f = 0 (corresponding to the edge of the j et) . In this case 
p'f is definitely not small, from Equation (D.18). This point shall be 
discussed further when the results for the air-water jet are considered. 

The Boussinesq approximation for scalar quantities is (47) 


(|)»U' 




(D.19) 


4> 'v' 


^t 9(}) 
a , 3r 


(D.20) 


which can be applied to model terms (c) and (d) Equation (D.14), with 
0 = f. Equations (D.3) and (D.4) can be used to model terms (i) and 

(j). 


Employing the results of the previous discussion. Equation (D.14) 
can be rewritten as follows: 


(a) 


pu 


9f 

9x 


(b) 


+ pv 


3f 

8r 


9x 


(c). 


(P 


(i) 

9p 

9x 

a 


9f 

9x 


(j) 

_9p ^ 

9r 9r 


9f , 1 

0 + — 


9x 


9r 


Cd) 
(rp ■ 


a, 


9£ 

9r 


(D.21) 


Neglecting higher-order terms and employing Equation (D.8), Equation 
(D.21) becomes 


ax --0 
pu ^ + pv 


9r 


1 9 

r 9r 


Cr 


y 


t 


a 


f 


9f 

9r 


-) = 0 


(D.22) 


The turbulent kinetic energy equation requires modeling additional terms. 
Employing the argument given by Gosman, et al., (47), the following 
equation is obtained after completion of the order of magnitude estimation 
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— 3k 


+ pv° 


3r 


1 ^ 

r 9r 


t 3k . 
°k 


= y. 


, 3u . 2 
^ 3r '' 


pe 


CD. 23) 


In arriving at this equation, a number of density fluctuation terms are 
ignored and other terms are modeled according to the methods employed 
for uniform property flow. In fact, the only difference betwee n Equ ation 
(D.23) and the constant density form involves retention of the p'v' 
correlation term on the left hand side of the equation. The justification 
for this largely rests on the success of the method during earlier studies 
of variable density flows (47). 


Similar simplifications are empl oyed for the conservat ion e quations for 
e and g. As in the case of the p’f correlation, the p*g’ correlation 
(which could be modeled using the probability density function) has also 
been ignored. The final results are: 


pu 


3e 

3x 


+ 


pv° 


3e 

3r 


C 


e 


P 

2 



13^ ^t 3 e ^ 
r ^ 3F"^ 


^t 

El 


£ ^ 9u ^ 2 

k ^'3r ^ 


(D.24) 


pu 



pv° 


3r 


1 ^ 

r 3r 


(r 




®i " 


3f ^2 


3r 



(D.25) 


The constants, 
turbulent 
a, , a and a 
dow/ 8 


C C , C and C , as well as the values for the 
^1' ^2 ^1 ®2 Prandtl/Schmidt Numbers, 

are drawn from the values obtained for uniform property 


Given the local values of k and e, the model is completed by obtaining 
the turbulent viscosity from- the following equation (47) 


Wt = % Pk^/e 


(D.26) 


where 


is also a constant specified from uniform property flow resultTs. 
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APPENDIX E 

EVALUATION OF PROPERTIES FOR DROPLET LIFE HISTORY MODEL 


E.l Thermal Conductivity 

The thermal conductivity of the gaseous mixture was calculated with the 
Mason and Saxena formulation of the Wassiljewa Equation (65). 


N I — X . X . “1 

» ■ . 


I 


j=l 




where 


- 1/2 

= 8 (1 - M 


”i . 1/2 rr , /i w Mi .1/4] 

-M-) L' " ^ V ^ -J 


(E.l) 


(E.2) 


Variation of the thermal conductivity of the gases was taken to be a 
linear function of temperature in the range of the calculations. 

E.2 Viscosity 

The viscosity calculations employed the method of Wilke as described by 
Reference (65) . 

For a mixture of N components 


U 


N 


I 

i=l 



(E.3) 


where (j) is given by Equation (E.2), 
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E.3 Gas Properties 


Viscosity 



Temperature (K) 

y (Ns/m^) 

Source 


Freon-11 

258.15 

.087 X 10”^ 

66 


303.15 

1.08 X 10"^ 

66 

Air 

240 

1.547 X 10"^ 

63 


300 

1.846 X 10 ^ 

63 


Thermal 

Conductivity 




Temperature (K) 

A (W/mK) 

Source 

Freon-11 

258.15 

5.892 X 10“^ 

66 


303.15 

7.799 X 10"^ 

66 

Air 

240 

2.145 X 10“^ 

63 


2.624 X 10 


300 


63 
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APPENDIX F 


EXPERIMENTAL DATA 


F. 1 Axial Variation of Quantities 
Isothermal single-phase jet 




x/d 

C O 


50 

.142 

108 

.065 

170 

.045 

340 

.022 

510 

.016 


Variable Density Single-Phase Jet 



x/ d 

u /u 
c o 

x/d 



50 

.351 

110 

.112 

108 

.163 

170 

.073 

170 

.096 

340 

.033 

340 

.058 

510 

.019 

510 

.041 




= 1 
o 
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Evaporating Spray 



x/d 

u /u 
c o 

x/d 

f /f ^ 
c o 

x/d 

(T -T )/(T -T ) 

' Q oo' ' O 


50 

.372 

75 

.268 

300 

.629 

78 

.275 

170 

.144 

340 

.465 

170 

.157 

340 

.076 

375 

.384 

267 

.108 

510 

.047 

425 

.337 

3 A 0 

.072 



450 

,294 

435 

.059 



510 

.262 

510 

.050 
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F.2 Radial Variation of Quantities 


Isothermal Single-Phase Jet 



x/d 

170 

340 


510 


r/x 

u/u 
' c 

r/x 

u/u^ 

r/x 

a/uj. 


026 

.948 

.019 

.972 

.021 

.936 

039 

.815 

.034 

.896 

.042 

.835 

055 

.783 

.049 

.835 

.064 

.754 

072 

.632 

.061 

.746 

.085 

.561 

088 

.576 

.078 

.637 

.097 

.400 

104 

.470 

.094 

.496 

.146 

.219 

118 

.364 

.109 

.418 

.169 

.123 

136 

.282 

.125 

.315 

.180 

.068 

164 

.141 

.141 

.238 

.209 

.025 

178 

.085 

.154 

.184 



199 

.042 

.174 

.116 





.185 

.062 





.203 

.034 
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Isothermal Single-Phase Jet 


x/d 170 

340 

510 

r / X 

k/uc 

r/ac 

k/u^ 

c 

r/x 

k/u^ 

c 

.015 

.0768 

.022 

.0802 

.023 

.0765 

.043 

.0778 

.042 

.0802 

.044 

.0765 

.062 

.0730 

.062 

.0771 

.066 

.0756 

.084 

.0561 

.083 

.0672 

.087 

.0642 

.104 

.0486 

.103 

.0568 

.105 

.0533 

.124 

.0304 

.125 

.0427 

.128 

.0383 

.146 

.0197 

.144 

.0321 

.147 

.0272 

.166 

.0124 

.166 

.0238 

.170 

.0190 

.183 

.0078 

.187 

.0155 

.192 

.0105 


Isothermal Single-Phase Jet 


x/d 170 

340 

510 

r/x 


r/x 

o 

r/x 


u’vVu*- 

c 

u’v’ /u 

c 

u'v’ /u"^ 
c 

.021 

.0038 

.024 

.0094 

.021 

.0062 

.033 

.0117 

.039 

.0142 

.042 

.0139 

.049 

.0158 

.058 

.0171 

.063 

.0189 

.063 

.0183 

.072 

.0182 

.085 

.0192 

.078 

.0188 

.084 

,0183 

.104 

.0175 

.094 

.0181 

.102 

.0172 

.125 

P140 

.109 

.0160 

.118 

.0150 

.146 

.0088 

.124 

.0134 

.135 

.0109 

.168 

.0053 

.138 

.0101 

.150 

.0082 

.186 

.0033 

.152 

.0069 

.179 

.0039 

.209 

.0015 

.170 

.0044 

.196 

.0022 



.186 

.0026 

.211 

.0009 



.207 

.0011 
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Variable Density Single-Phase Jet 


x/d 170'^ 340 510 


r/x 

u/ u 

c 

r/x 

u/u 

c 

r/ X 

u/u 

c 

.019 

.955 

.014 

.988 

.022 

.971 

.036 

.852 

.028 

.945 

.044 

.867 

.054 

.759 

.059 

.762 

.063 

.749 

.072 

.621 

.076 

.664 

.084 

.608 

.090 

.525 

.090 

.550 

.105 

.474 

. 108 

.378 

.105 

.448 

. 124 

.344 

.126 

.286 

.121 

.357 

.145 

.242 

.145 

.193 

.133 

.274 

.166 

.176 

.163 

.130 

.151 

.195 

.186 

.098 

.182 

.082 

.167 

.146 

.207 

.056 

.199 

.045 

.184 

.087 

.229 

.043 

.217 

.023 

.198 

.066 
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Variable Single-Phase Density Jet 


x/ d 

170 

340 


510 


r/x 


r/x 


r/x 


.037 

.880 

.018 

.991 

.023 

.979 

.054 

.782 

.028 

.955 

.046 

.845 

.074 

.638 

.042 

.983 

.067 

.768 

.093 

.536 

.062 

.782 

.085 

.627 

.107 

.412 

.074 

.658 

.106 

.491 

.126 

.310 

.092 

.565 

.126 

.365 

.145 

.233 

.106 

.462 

.146 

.283 

.150 

.163 

.119 

.344 

.165 

.136 

.182 

.099 

.134 

.290 

.190 

.121 

.200 

.067 

.151 

.221 

.210 

.087 



.168 

.170 





.184 

.119 





.199 

.094 
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Variable Density Single-Phase Jet 


x/ d 

170 

340 


510' 


r/ X 

k/u^ 

c 

r/x 

k/u^ 

c 

r/ X 

V/-2 

k/u 

c 

0 

.0726 

0 

.0789 

0 

.0813 

.021 

.0728 

.022 

.0816 

.025 

.0784 

.045 

.0735 

.044 

.0810 

.049 

.0783 

.062 

.0663 

.062 

.0788 

.068 

.0765 

.082 

.0592 

.084 

.0656 

.091 

.0647 

.106 

.0487 

.105 

.0553 

.111 

.0530 

.123 

.0359 

.125 

.0451 

.130 

.0414 

.140 

.0232 

.145 

.0306 

.150 

.0273 

.166 

.0104 

.167 

.0205 

.170 

.0165 

.187 

.0048 

.187 

.0124 

.192 

.0048 
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Variable Density Single-Phase Jet 


x/d 170 340 510 


r/x 

u' v’ /u 

c 

r/x 

1 ? /“2 
U V /u 

c 

r/x 

u’v'/u 

.019 

.0076 

.018 

.0067 

.024 

.0091 

.036 

.0119 

.029 

.0108 

.042 

.0152 

.053 

.0122 

.060 

.0183 

.062 

.0188 

.072 

.0127 

.076 

.0186 

.084 

.0191 

.092 

.0121 

.088 

.0175 

.104 

.0172 

.108 

.0102 

.121 

.0148 

.124 

.0128 

.124 

.0078 

.132 

.0107 

.145 

.0093 

.145 

.0051 

.152 

.0076 

.165 

.0062 

.167 

.0031 

.168 

.0049 

.186 

.0025 

.183 

.0016 

.185 

.0031 

.228 

.0009 



.200 

.0022 
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Evaporating Spray 


x/d 

170 

340 


510 


r/ X 

u/ u 

c 

r/x 

u/u 

c 

r/x 

u/u 

c 

.033 

.905 

.016 

.947 

.022 

.971 

.049 

.817 

.028 

.913 

.043 

.878 

.063 

.693 

.045 

.829 

.065 

.749 

.079 

.788 

.064 

, 662 

.084 

.585 

.096 

.476 

.077 

.608 

.094 

.450 

.110 

.355 

.092 

.496 

.123 

. 328 

.126 

.280 

.107 

.386 

.145 

.210 

.156 

.137 

.120 

.283 

.166 

.113 

.188 

.039 

.137 

.207 

.188 

.061 



.152 

.160 





.163 

.092 
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Evaporating Spray 


x/ d 

170 

340 


510 


r/x 


r/x 


r/x 


.025 

.981 

.023 

.957 

.016 

.992 

.049 

.876 

.043 

.915 

.034 

.939 

.076 

.723 

.063 

.788 

.048 

.921 

.095 

.568 

.085 

. 654 

.062 

.817 

.117 

.452 

.107 

.523 

.078 

.705 

.132 

.323 

.124 

.455 

.094 

.618 

.166 

.188 

.162 

.247 

.125 

.428 

.189 

.102 

.186 

.138 

.141 

.364 


.155 .248 


172 


187 
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Evaporating Spray 


x/d 

170 

340 


510 


r/x 

k/u^ 

c 

r/x 

k/u^ 

c 

r/x 

k/u^ 

c 

0 

.0701 

0 

.0743 

0 

.0791 

.022 

.0708 

.020 

.0750 

.024 

.0780 

.041 

.0685 

.041 

.0752 

,043 

.0777 

.061 

.0671 

.061 

.0748 

.065 

.0736 

.082 

.0613 

.082 

.0652 

.087 

.0662 

.102 

.0517 

.102 

.0548 

.105 

.0563 

.122 

.0385 

.122 

.0422 

.126 

.0446 

.141 

.0255 

.140 

.0293 

.146 

.0318 

.163 

.0137 

.163 

.0175 

.166 

.0207 

.182 

.0026 

.182 

.0027 

.185 

.0108 
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Evaporating Spray 


x/d 170 340 510 

r/x u’v’/u^ r/x u’v'/u^ r/x u'v'/u^ 

c c c 


.017 

.0016 

.021 

.0061 

.017 

.0067 

.035 

.0057 

.033 

.0122 

.032 

.0135 

.048 

.0088 

.041 

.0153 

.048 

.0176 

.063 

.0112 

.064 

.0144 

.062 

.0193 

.078 

.0132 

.085 

.0113 

.078 

.0190 

.094 

.0135 

.104 

.0102 

.108 

.0141 

.109 

.0120 

.126 

.0061 

.112 

.0080 

.127 

.0089 

.147 

.0039 

.142 

.0047 

.142 

.0063 

.168 

.0025 

.158 

.0022 

.157 

.0042 

.187 

.0012 

.173 

.0013 


.174 .0024 

.189 .0020 
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Evaporating Spray 


x/d ~34Q~~ ’ ^ 


r/x 

(T-T )/(T -T ) 

oo' ' Q oo' 

r/x 

(T-T )/(T^-Tj 

OO Q W 

.023 

.907 

.021 

.931 

.047 

.698 

.042 

.783 

.075 

.543 

.062 

.663 

.094 

,387 

.082 

.510 

.116 

.271 

.104 

.426 

.141 

.223 

.125 

.318 

.163 

.128 

.145 

.198 

.188 

.063 

.167 

.125 



.186 

.085 
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